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To my mother and father 




Preface 



The interaction of an electron beam with a solid target has been studied since 
the early part of the past century. Since 1960, the electron-solid interaction 
has become the subject of a number of investigators’ work owing to its funda- 
mental role in scanning electron microscopy, in electron-probe microanalysis, 
in Auger electron spectroscopy, in electron-beam lithography and in radiation 
damage. The interaction of an electron beam with a solid target has often 
been investigated theoretically by using the Monte Carlo method, a numeri- 
cal procedure involving random numbers that is able to solve mathematical 
problems. This method is very useful for the study of electron penetration 
in matter. The probabilistic laws of the interaction of an individual electron 
with the atoms constituting the target are well known. Consequently, it is 
possible to compute the macroscopic characteristics of interaction processes 
by simulating a large number of real trajectories, and then averaging them. 
The aim of this book is to study the probabilistic laws of the interaction of 
individual electrons with atoms (elastic and inelastic cross-sections); to in- 
vestigate selected aspects of electron interaction with matter (backscattering 
coefficients for bulk targets, absorption, backscattering and transmission for 
both supported and unsupported thin films, implantation profiles, secondary- 
electron emission, and so on); and to introduce the Monte Carlo method and 
its applications to compute the macroscopic characteristics of the interac- 
tion processes mentioned above. The book compares theory, computational 
simulations and experimental data in order to offer a more global vision. 

I wish to thank Antonio Miotello (Trento University) and Francesc Salvat 
(Barcelona University) for their invaluable and stimulating comments. I am 
also very pleased to acknowledge Angela Lahee (Springer- Verlag) for her con- 
structive cooperation. I appreciate the conscientiousness of Douglas Meekison 
(Springer- Verlag) who read the entire manuscript and bettered it consider- 
ably. I am indebted to Hayley Stead for her skilful technical assistance in 
checking and improving the quality of the English of the book. Warm thanks 
are due to my wife and children for their immense patience and understanding 
during the period when this monograph was written. 



Povo, January 2003 



Maurizio Dapor 
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1 Introduction 



1.1 Motivation 

Studies of low-energy backscattered electrons, secondary electrons and ab- 
sorbed electron depth distribution, have all been of great interest for many 
analytical techniques. Indeed, a better understanding of the collision events 
occurring in the surface layers before the emission of backscattered and 
secondary electrons, should allow a more general comprehension of surface 
physics. The problem of backscattered and secondary electron emission from 
solids irradiated by a particle beam is of crucial importance, especially in 
connection with the analytical techniques that utilise secondary electrons to 
investigate chemical and compositional properties of solids in the near surface 
layers, i.e. Auger electron spectroscopy and X-photoelectron spectroscopy. 
The electron backscattering coefficient is a quantity which, for slow primary 
electrons and for supported and unsupported thin films, requires theoretical, 
numerical and experimental investigation. The energy spectra of secondary 
electrons emitted by solid targets and stimulated by particle beams are very 
complicated because many features appear in the spectra due to the different 
collision processes involved prior to low-energy secondary electron emission. 
Accurate knowledge of the depth distribution of trapped electrons in dielec- 
tric materials (biological materials, ceramics, glass) is necessary to solve the 
diffusion equation in order to investigate the charging phenomena, while the 
study of positron depth distribution in solids is required for positron annihi- 
lation spectroscopy, a technique that allows non-destructive investigations of 
structural defects at surfaces and interfaces. 

As a consequence, the backscattered electrons and positrons, the sec- 
ondary electrons, and the depth distribution of the electrons and positrons 
absorbed by solid targets, all need to be theoretically and numerically inves- 
tigated. 

These, and other related quantities, can be accurately calculated using 
Monte Carlo simulations. 

This book aims to provide the reader with the physical and mathematical 
instruments necessary to realise a Monte Carlo code in order to simulate the 
quantities quoted. But before proceeding, we would like to briefly add some 
words about the physical quantities just introduced. 
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1 Introduction 



1.1.1 Backscattered Particles 

When an electron or a positron beam impinges on a solid, some particles of 
the beam scatter back and emerge from the target surface. These particles 
are known as backscattered electrons. Because of their penetration below the 
surface and the resulting loss of small amounts of energy through ionizations, 
electron excitations and plasmon emissions, the backscattered electrons are 
not reflected without dissipation of energy. A fraction of their energy is indeed 
lost into the solid before emerging. 

When the target is a surface film (i.e. a thin film deposited on a substrate 
of a different material), the number of backscattered electrons depends on 
the film thickness, on the electron primary energy, and on the materials con- 
stituting the film and the substrate. 

In simulating the processes of penetration and backscattering of electrons 
impinging on supported thin films with the Monte Carlo method, one has to 
take into account not only the interface between the film and the substrate 
but also the change in the scattering probabilities per unit length in passing 
from the film to the substrate and vice versa. 

The fractions of electrons backscattered from unsupported thin films, solid 
targets, and surface films of various thickness and deposited on different sub- 
strates, can be accurately calculated using Monte Carlo simulations. 

1.1.2 Secondary Electrons 

The concepts which form the basis for a theoretical understanding of the 
problem of secondary electron emission have been analysed in Wolff’s theory 
[1]. However, such a theory must include some simplification which may be 
valid in only a limited number of situations in order to attain analytical 
results. 

To overcome these limitations, we have therefore set out, in this book, 
to describe a Monte Carlo procedure to accurately calculate the number of 
secondary electrons emitted from a solid irradiated with positron or electron 
beams. 

1.1.3 Depth Distribution of Trapped Particles 

When an insulator is subject to electron irradiation, the injected electrons 
cannot be definitely trapped but must instead recombine with positive 
charges left near the irradiated surface when secondary electrons are emitted: 
this is justified on the basis that dielectric breakdown is not observed dur- 
ing specific experiments of electron-irradiation of insulators. The dynamics 
of the absorbed electrons depend on a number of parameters: the fraction of 
trapped electrons, the space-charge distribution, the mobility, and the num- 
ber of secondary electrons emitted from the region near the surface of the 
dielectric. The time evolution of the surface electric field can be studied by 
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integrating the continuity equation for the relevant transport processes of the 
injected charge by adopting, as the charge source term, the distribution of 
the absorbed electrons as obtained by a Monte Carlo simulation. 

The study of the interaction of positron beams with solid targets has been 
approached by various researchers, due to its importance for positron anni- 
hilation spectroscopy. This technique allows non-destructive investigations of 
the structural defects of surfaces and interfaces: in particular, information is 
provided about the nature and distribution of point defects in solid materials. 
The solution of the diffusion equation (necessary to obtain the fractions of 
incident positrons annihilated at different depths within the target) requires 
knowledge of the positron stopping profile, i.e., the initial depth distribution 
of the thermalized positrons that can be calculated through Monte Carlo sim- 
ulation. We also note that the transmission of positrons is of great interest 
because, once backscattering is known, it allows calculation of the total frac- 
tion of particles absorbed by the target as a function of depth and primary 
energy. The calculation of the transmission of particles through thin films is 
also a typical application of Monte Carlo simulations and will be described 
in the last chapter. 



1.2 Prom the Dirac Equation 

to the Monte Carlo Simulation 

In the previous section, several examples of physical problems concerning 
the interaction of electrons and/or positrons with solids have been given. We 
have also noted that these problems are often solved by using Monte Carlo 
simulation, a computational tool based on the use of random variables. 

An accurate knowledge of the elastic and inelastic cross-sections is needed 
to study the electron and positron penetration in solids. 

The calculation of the scattering cross-sections, necessary for the Monte 
Carlo simulations, requires the study of the electron spin and the Dirac equa- 
tion. This important equation will be described in Chap. 2. The following 
chapters will be dedicated to the problem of calculating cross-sections and 
to the description of the physical quantities involved in the electron-solid 
interaction that we are interested in investigating. The last chapter describes 
the Monte Carlo codes. 

Numerical procedures to calculate the differential elastic scattering cross- 
sections necessary in the Monte Carlo simulations are described in Chap. 3. 

The codes for the calculation of the elastic scattering cross-section are 
based on the quantum-relativistic partial- wave expansion method. The phase 
shifts (necessary to compute the differential elastic scattering cross-section) 
can be calculated by numerically solving the Dirac equation for a central 
electrostatic field up to a large radius in which the atomic potential is neg- 
ligible. The atomic potential used here is that of Hartree-Fock for atomic 
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1 Introduction 



numbers smaller than or equal to 18, and that of Dirac-Hartree-Fock-Slater 
for atomic numbers greater than 18. 

The inelastic scattering can be quantified by using the inelastic mean free 
path and the stopping power, a quantity related to the probability of energy 
loss per unit distance travelled by the electron within the solid. An electron 
can lose a large fraction of its energy in a single collision: nevertheless the 
continuous slowing down approximation is generally accepted. In this approx- 
imation, the electron is assumed to continuously dissipate its energy whilst 
travelling within the solid. In this book, the inelastic mean free path and 
the stopping power are described following the method proposed by Ashley 
[2, 3]. Numerical procedures to calculate the inelastic mean free path and the 
stopping power, necessary in the Monte Carlo simulations, are described in 
Chap. 4. 

Selected aspects of electron interaction with matter (backscattering co- 
efficients for bulk targets, absorption, backscattering and transmission for 
both supported and unsupported thin films, implantation profiles, secondary- 
electron emission, and so on) are described in Chap. 5. 

The Monte Carlo codes described in Chap. 6 give results which agree 
well with the experimental data in the electron and positron energy range 
from 1 000 to 30 000 eV. The numerical and theoretical results presented in 
this book have been utilised by several authors in their simulations (see, for 
example, [4, 5, 6]). 

The mathematical and physical requirement are kept within the topics 
covered by a calculus course and by an introductory quantum mechanics 
course, respectively. The appendix provides some useful tools, i.e. a summary 
of the theory of matrices and operators, the Dirac notation (ket and bra 
vectors) and a description of the special functions of the mathematical physics 
(Legendre polynomials, associated Legendre functions, Bessel functions and 
spherical harmonics) . 
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2 The Spin of the Electron 



2.1 The Spectrum of Angular Momentum 

Let us introduce the angular momentum operator in the treatment of single 
particle quantum systems [1, 2, 3, 4]. Let us indicate the electron mass by 
m, the electron position by r, the electron energy by E and the electron 
momentum by p. In quantum mechanics we assume the following correspon- 
dence rules relating the differential operators (Appendix A) and the physical 
quantities: 



'’'Iff 


(2.1) 


-ifiV . 


(2.2) 



Here Ti = h/2TT and h = 4.136 x 10“^^ eV sec is the Planck constant. The 
differential operators act on wave functions that are square-integrable com- 
plex functions in a Hilbert space. Now, if we consider the components of the 
electron orbital angular momentum L = r x p, using the definition of L it is 



possible to see that 




— iflLz ; 


(2.3) 


\Ly^Lz\ — iTiLx ; 


(2.4) 


\Lzi L/x\ — iTiLy . 


(2.5) 


The commutator is defined for any pair of linear operators A and B 


as 


[x4, B] = AB - BA . 


(2.6) 


In order to introduce the intrinsic angular momentum, we can begin by 
generalising the properties of L, saying that a linear operator J is an angular 
momentum (orbital or intrinsic) if 


[Jxt Jy\ iTlJz ; 


(2.7) 


[t/y, J z\ iTtJx t 


(2.8) 


\Jz-> Jx\ — ■ 


(2.9) 
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2 The Spin of the Electron 



As the reader may easily verify, the two operators and Jz commute, i.e. 

[J.,J'] = 0, (2.10) 

so they possess at least one basis of eigenvectors in common. As a conse- 
quence, the corresponding physical quantities may be measured together and 
simultaneously with an arbitrary precision. Using the Dirac notation (Ap- 
pendix B) and using \jm) to indicate a basis of orthonormal eigenvectors 



that and Jz have in common, the eigenvalues j{j + 1) and hm corre- 
spond respectively to and Jz, whereby 

= n^j{j + l)\jm) , (2.11) 

Jz\jm) = . (2-12) 

In order to proceed, let us now define the two operators J_ and J+ by 

J± = Jx ^ iJy , (2-13) 

which have the following properties: 

4 = , (2.14) 

[Jz,J±] = ±hJ± , (2.15) 

J^J— = J^ Jy flJz ) ( 2 - 16 ) 

J-J+ = Jl + j'^-hJz , (2.17) 

[J+,J-] = 2hJz , (2.18) 

{J+,J_} = 2(j2 + J2) . (2.19) 

Here we have introduced the anticommutator, defined for any pair of linear 
operators A and B by 

{A,B} = AB + BA. (2.20) 

As a consequence, 

= Jl+\{J+J- + J-J+) ■ ( 2 . 21 ) 



Let us now calculate the norms of the vectors J-\jm) and J+|jm), i.e. 
(jm| 4 J_ |jm) = {jm\J+J-\jm) = -I- 1) — m{m — 1] > 0 , 

( 2 . 22 ) 

{jm\j\_J+\jm) = {jm\J-J+\jm) = lf[j{j + 1) - m{m + l] > 0 . 

(2.23) 



As the norms are positive, 
- j <m< j . 



(2.24) 
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Note that m = —j if and only if J_|jm) = 0, while m = j if and only if 
= 0. As commutes with J±, then 

J^J±\jm) = J±J^\jm) = h^j{j + l)J±\jm) . (2.25) 

From (2.15), we obtain 

JzJ±\j'm) = J±Jz\j'm) ± hJ±\jm) = h(m± l)J±\jrn) . (2.26) 

The last two equations tell us that the vectors J± \jm) are common eigenvec- 
tors of and Jj. On the other hand, comparing the last equation with 

Jz\jm±l) = h{m±l)\jm±l) , (2.27) 

we can see that the eigenvectors J±\jm) must be proportional to \jm ± 1). 
Therefore, taking into account the values of their norms that we have just 
calculated, we can write that 

J± \jm) = Tl\/j{j + 1) - m(m± 1)1 jm ± 1) . (2.28) 



Following the last equations, simple considerations about the integer number 
of steps p necessary to go from m = —j to m = j using the operator J+ 
i~j +P = j) allow us to conclude that j must be an integral or half-integral 
non-negative number. 





,00 , 



(2.29) 



and that the only possibles values of m are integral and half-integral numbers 



m = 0,±-,±l,±-,±2,±-, . . 



, ±oo . 



(2.30) 



In conclusion, if hj{j + 1) and hm are, respectively, the eigenvalues of and 
Jz then j must be an integral or half-integral non-negative quantity, and the 
only possible values of m for any given j are the 2j -|- 1 numbers 



+ ■■■ ■ 



(2.31) 



2.2 The Spin of the Electron 

As we are interested in electrons and in their intrinsic angular momentum, 
or spin, we shall now focus our attention on the case j = 1/2. Let us simplify 
the equations by introducing the following notation: 

I-) = |l/2 - 1/2) , (2.32) 

1+) = |l/2 + 1/2) . (2.33) 

Let us indicate the spin operator by S and introduce the two-dimensional 
spin eigenspace defined by the two eigenvectors |— ) and |-|-), possessing the 
following properties: 
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(+l-) = H+) = o, 


(2.34) 


(-!-> = (+l+) = i, 


(2.35) 




(2.36) 




(2.37) 


50 

T 

II 

1 

to| St 


(2.38) 


'S'z +) “ 2 ■ 


(2.39) 



The general case of a spin- 1/2 state ^ is a linear superimposition of the basic 
eigenvectors |— ) and |-|-), whereby 

^ = A\+) + B\~), (2.40) 

where the coefficients A and B are complex numbers. \A\^ is the probability 
of finding the electron in the state of “spin up” along the z axis; \B\’^ is the 
probability of finding the electron in the “spin down” state. It follows that, 
from the condition of normalisation of 



|Ap-h|Bp = l. (2.41) 

Using the equations 

^-|-) = 0, (2.42) 

S_\+) = n\-), (2.43) 

5+I-) = ;i|+) , (2.44) 

5+1+) = 0 , (2.45) 



we can calculate the matrix elements of 5+ and 5_, as follows: 

( 0 J ) ’ ( 2 - 46 ) 

5_ = J 0 ) • (2.47) 

Let us now define the Pauli matrices ay and cr^ by 

5 = |<T. (2.48) 

The representation of the Pauli matrices in the basis {|— ), |+)} can be easily 
obtained, observing that 
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S, = ^{S. + S+), 

Sy = - S+) . 

The Pauli matrices in the basis {|— ), |+)} are 




(2.49) 

(2.50) 

(2.51) 

(2.52) 



) . (2.53) 

Using this particular representation or their definition, we can see that the 
Pauli matrices are Hermitian and have the following properties: 

al = al = a^, = l, (2.54) 

= i , (2.55) 

O'xO'y (Ty(7x — z ; (2.56) 

CTzCTx - <7x<7z = , (2-57) 

ayaz - <TzC7y = 2iax , (2.58) 

(7x<Jy + CTyax = 0 , (2.59) 

CTzCTx + O-xCTz = 0 , (2.60) 

o-yCTz + azCTy = 0 , (2-61) 

tr Ua; = tr CTy = tr CTz = 0 , (2.62) 

det ax = det ay = det cr^ = — 1 . (2.63) 
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When applied to the energy of the particle they give the non-relativistic 
Hamiltonian operator of a particle in a field, 



^ = -^V2 + H(r) , 



(2.67) 



and the Schrodinger equation describing the time evolution of the wave func- 
tion, 






2m 



■V{r) 






( 2 . 68 ) 



or 



at 



(2.69) 



The Dirac equation is the appropriate equation to treat the electron spin 
and relativistic effects. The square of the energy if of a relativistic free elec- 
tron is given by 

, (2.70) 



where c is the speed of light, m is the rest mass of the electron and p is the 
electron momentum. In the following, all equations will be written with a 
choice of units such that 



Ti = c = 1 , 



(2.71) 



in which time has the dimensions of length, the electron charge e is a dimen- 
sionless quantity such that 



e 



2 



1 

m ’ 



(2.72) 



and the square of the energy of a relativistic free electron becomes 



E^ = p^ 



(2.73) 



Note that through homogeneity considerations, we can easily re-establish the 
general equations. 

For the free electron, the Dirac Hamiltonian is given by 



E[ = a. ■ p + (3m = —ia ■ V -f (3m , (2.74) 

in which a = {ax,ocy, ocz) = (a^, a^, a^) and (3 are four Hermitian operators. 
Therefore, 

d 

t— (() = [—ia. ■ V -I- (3m\ (p , (2-75) 

where </> is the wave function. The square of the Dirac Hamiltonian must 
obviously be equal to p^ + m^-. 

{axPx + ctyPy + azPz T (3mf = p^ + m^. 



(2.76) 




2.4 The Solution for Free Particles 
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and therefore 

+ y^(a-^Q:^ + a’^a^)p’p^ + P + Pa^)p^ 

j j<k 3 

= ttP + '^p>^ .(2.77) 

3 

As a consequence, the four operators and P must satisfy the following 
relations: 



a^'2 =p^ = l^ 

= 0 for j fc , 

a^P + Pa^ = 0 . (2.78) 



As they must obey these relations, the operators ai and P obviously can- 
not be numbers; instead they must be Hermitian matrices operating on 
a multi-component column vector p (which we call a spinor). The lowest- 
dimensionality matrices that obey the relations (2.78) are 4x4. If 
((Tx, cTy, az) = (cr^, cr^, <J^) are the 2x2 Pauli matrices, and I the 2x2 identity 
matrix, then the following 4x4 matrices satisfy the conditions expressed by 
(2.78): 



a 



3 



0 \ 

0 j ’ 



(2.79) 



P = 



I 

0 




(2.80) 



In the presence of an electromagnetic field described by the four-potential 
(v?, A), the Dirac Hamiltonian becomes {E ^ E — e<p, p ^ p — eA) 



H = eifi + a - {p — eA) + Pm 



(2.81) 



so that the Dirac equation describing a spin-1/2 particle in the presence of 
an electromagnetic field is 

d \ 

— — epj — a • (— iV — eA) — /3m p = 0 , (2.82) 

where p{r,t) is a four-component spinor depending on r and on t. 



2.4 The Solution for Free Particles 

In the case of free electrons, the Hamiltonian does not depend on r and t, so 
that we can express p as & plane wave: 

p = u{p) exp[i(p • r — Et)\ . 



(2.83) 
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As (j) must obey the free Dirac equation, then the function u{p) must be a 
four-component spinor independent of r that satisfies the equation 

{ot ■ p + Pm)u{p) = Eu{p) . (2.84) 

This is equivalent to 

{ {E-m)ui -PzU 3 -{px-ipyjui =0 

{E - m)u 2 ~{Px + iPy)u3 +PzU4 =0 

-PzUl -{Px-ipy)u 2 +{E + m)u3 =0 

-{Px + ipy)ui +PzU2 +{E + m)u 4 =0 

or, by choosing the 2 : axis in the direction of p {px = Py = 0, Pz = p), 



{E — m)ui 


-PU3 


= 0 




1 

to 


+PU4 = 0 


-pui 


+{E + m)u3 


= 0 




+PU2 


+ {E + m)u4 = 0 



In order to find the eigenvalues of the energy, we have to solve the following 
equation: 



E — m 0 

0 E — m 

—p 0 

0 p 



—p 0 

0 p 

E + m 0 

0 E + m 



= 0 



Upon solving this last equation, we can deduce the following doubly degen- 
erate eigenvalues for the energy of the free Dirac particle: 



E± = ±\/p2_^m2 . (2.85) 

Note that the spin operator in the direction of p is the operator CTz/2 = 
—iaxayj2, which commutes with the Hamiltonian iJ = azP+jSm, so that the 
solutions we are seeking correspond to the eigenvectors which are common to 
E[ and (Tz/ 2. In order to find the first eigenvector, let us just impose rti = 1 
and U 2 = 0 and normalise to unity {u^u = 1). In such a way, we find that 



U^E+{p) 



1 


/ 


1 


\ 


1 -{-m 




0 




/ 2if+ 




p/(E+ + m) 






V 


0 


/ 



( 2 . 86 ) 



where e+ (p) is the eigensolution corresponding to a particle with spin up 
and positive energy (i.e. a spin-up electron). The other eigenvectors can be 
obtained by proceeding in a similar way. The solution corresponding to a 
spin-down electron is 



1 


( ^ \ 


E+ + m 


1 


/ 2A+ 


0 




\ -p/{E+ + m) J 



UlE+{p) 



(2.87) 
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while the two spinors corresponding to the positron (negative energy) are 



U^E^ip) 



UlE.{p) 



I 


/ -p/{E+ + m) \ 


E+ + m 


0 


/ 2if+ 


1 




1 0 / 



1 




/E+ + m 


p/{E+ + m) 


/ 2if+ 


0 




1 1 / 



( 2 . 88 ) 



(2.89) 



Note that in the non-relativistic limit the solutions corresponding to the 
electron (positive energy) have the two components ui and U 2 much greater 
than t63 and so that, for w <C c, the four-component Dirac spinors reduce 
to the two-component spinors of the Pauli theory. 



2.5 The Dirac Equation in a Central Potential 

In order to appropriately treat the quantum-relativistic scattering theory, we 
need to know the form assumed by the Dirac equation for an electron in 
the presence of a central electrostatic field described by a central potential 
ep{r) = V{r). 

Let us first introduce the operator /C defined by 



K. = + L) , (2.90) 

where L is the electron orbital angular momentum. For an electron in a 
central electrostatic field, it is possible to show that 



^dL dcr 


(2.91) 


dt dt 


As a consequence, the total angular momentum, defined as J = 
is a constant of the motion. On the other hand. 


L+{l/2)a, 


J^-L^ = a-L+^, 


(2.92) 


so we can conclude that 




/C = /3(J2-L2 + i/4) , 


(2.93) 



1C commutes with H and is, as a consequence, a constant of the motion. 
Let us now define the radial- momentum operator pr, where 

_ .19 r ■ p — i 



r 



(2.94) 
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and introduce the radial component of the a operator, where 

ct • V 

ar = ; (2.95) 

r 

this obeys the relation 

al = l. (2.96) 

For any pair of vectors a and b, the following equations hold: 

{(T ■ a){cr ■ b) = a ■ b + i(T ■ a X b (2.97) 

and 



(a • a){a ■ b) = {cr ■ a)((T ■ b) . 
As a consequence, 

(a ■ r){a ■ p) = rpr + i(3JC . 
This last equation is equivalent to 

if}JC \ 



{OL- p) = ar[pr + 



r y 



As a result, the Dirac equation with the Hamiltonian 
ifiJC \ 



H = ar \Pr + 



r ) 



+ /3to + V (r) 



can be rewritten as the following: 

i!31C \ 



Oir [Pr + 



r J 



+ (3m + V{r) 



(j) = E((> 



(2.98) 

(2.99) 

( 2 . 100 ) 

( 2 . 101 ) 

(2.102) 



The operators (3, /C, and Jz are mutually commuting. In the following, X 
indicates an eigenvector common to those operators, so that 



PX = X, (2.103) 

K.X = -kX , (2.104) 

L'^X = l{l+l)X , (2.105) 

J,X = mjX. (2.106) 



Introducing the function y, where 

y = -arX, (2.107) 

we can observe that it has the following properties: 

X = -ary, (2.108) 

ICy = -ky, (2.109) 

(3y = ar(3X = -y . (2.110) 
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Let us now consider the following linear combination of X and y: 



(j)= F{r)y + iG{r)X , 



( 2 . 111 ) 



which is an eigenvector common to H, K, and Jz and thus the spinor we are 
looking for. 

Our objective is to determine the functions F{r) and G{r). The eigenval- 
ues of are (j±l/2)(j±l/2-|-l). As /C = /3(J^ — T^-l-1/4), the eigenvalues 
of 1C for the case j = 1 + 1/2 (spin up) are therefore given by 

k = — + -^ = —{I + 1) . (2.112) 

In the other case, where j = / — 1/2 (spin down), the eigenvalues of 1C are 

fc=(j + 0=/- (2.113) 



We are now able to find the equations corresponding to the radial be- 
haviour of the functions F and G. In order to do this, let us consider the 
Dirac equation (2.102) and observe the following: 



arPrF{r)y = i 



dF{r) 

dr 




(2.114) 



iarPrG{r)X 



dG{r) ^ G{r) 
dr r 



(2.115) 



tarfdlC ^ = --F{r)kX , (2.116) 

r r 

^-^C^iG{r)X = --G{r)ky , (2.117) 

r r 

/3mF{r)y = -mF{r)y , (2.118) 

(3miG{r)X = imG{r)X . (2.119) 



Hence, since X and y belong to different eigenvalues of f3 and are, conse- 
quently, linearly independent, the fundamental equations of the theory of the 
elastic scattering of electrons (and positrons) by atoms are 

[E + m- V{r)]F{r) + -h ^-^G(r) = 0 , (2.120) 

dr r 

-[E-m- F(r)]G(r) -h + ^-^F{r) = 0 . (2.121) 
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3 Elastic Scattering 



3.1 The First Born Approximation 



Let us begin our investigation by considering the non-relativistic problem of 
electron-atom elastic scattering, i.e. the scattering of a beam of electrons by 
a central potential V{r). 

The solid angle (W depends on the scattering angles [9, 9 + d9] and the 
azimuthal angles [</>, 4> + d(j)] . The differential elastic scattering cross-section 
da/df2 is defined as the ratio between the flux of particles per unit of time 
emerging in the solid angle df2 (divided by df2) and the incident flux. 

The flux of particles per unit time emerging after the collision in the 
solid angle dfi depends on the component jr of the current density in the 
outgoing direction from the centre of the atomic nucleus. The number of 
electrons emerging in the solid angle df2 per unit time is given by df2 
(note that dfi is the cross-sectional area normal to the radius). 

Let us consider a beam of incident electrons in the direction z normalised 
to one particle per unit of volume. Let K = mv/Ti be the electron momentum 
in the z direction, where v is the electron velocity, m the electron mass and 
h, the Planck constant divided by 2tt. This beam can be represented by the 
plane wave exp(iiLz). 

Since the incident beam has been normalised to one particle per unit of 
volume, then the electron velocity v is the incident flux. As a consequence. 






da jrT^ df2 



V dfi 



Jrr- 

V 



(3.1) 



At a large distance from the atomic nucleus, the potential V{r) is negli- 
gible and the scattered particles can be described by a spherical wave, i.e. a 
function which is proportional to exp(f ATr) /r. If /(0, (p) is the constant of pro- 
portionality (scattering amplitude), then the wave function W{r) of the whole 
scattering process (i.e. of the incident and the scattered electrons) satisfies 
the boundary conditions 

exp(f ATr) 



d'{r) 



exp(iATz) -k f{9,4>)- 



(3.2) 



The electron position probability density P is given by j'f'p = P*P, and the 
current density j{r,t) is 
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2m 



(3.3) 



Let us now calculate the radial component of the current density jr‘. 



Jr — 



no,4>) 



exp(iKr) d 



r dr 
exp(-iKr) d 






exp(— zitTr) 



dr 



no A) 



exp(iXr) 



v\no,ci)n^ 



(3.4) 



Comparing this equation with (3.1), we can see that the differential elastic 
scattering cross-section is the square of the modulus of f{9,(j))'- 



non) = ^ = \f{o,4>)? 



(3.5) 



In order to proceed, let us introduce the first Born approximation. 

The first Born approximation is a high-energy approximation. If E is the 
incident electron energy, e the electron charge, oq the Bohr radius and Z the 
target atomic number, the first Born approximation is quite accurate if 



E:$>—z^ 

Zao 



(3.6) 



Let us now introduce the Green function and the integral-equation approach. 
Starting from the Schrodinger equation. 



2tt7 

(V^ -h K‘^)E{r) = -^V{r)E{r) 
n 



(3.7) 



with the boundary condition expressed by (3.2), it is possible to show that 
this is a problem equivalent to the following integral equation: 



E{r) = exp{iKz) + J d^r' g{r,r')V{r')'E{r') , 

exp{iK\r — r'|) 



in which 

g{r,r') = - 



(3.8) 



(3.9) 



47r|r — r'\ 

is the Green function of the operator -I- K^. As is known, this operator 
satisfies the equation 

+ K‘^)g{r,r') = 5{') — r') , (3.10) 



where 5{r — r') is the Dirac delta function. 

Let us apply the operator -I- to the function E{r) defined by the 
integral equation (3.8): 

+ K'^)E{r) = + K'^)exp{iKz) 

OjY} C 

+ d^r'iV^ + K‘^)g{r, r')V{r')E{r') . (3.11) 



3.1 The First Born Approximation 
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The application of the operator to the plane wave exp(iiT 2 :) gives 
^2 

V^exp(iitTz) = exp(iitT2;) = —K^e-xp{iKz) 



and, as a consequence, we can write 
(V^ + K^) exp(iifz) = 0 . 
Therefore, 



(3.12) 

(3.13) 



(V^ + K-2)!f (r-) = 



[ d^r'{'V^ + K^)g{r,r')V{r')^{r') 
U J 

[ d^r'5{i — r')V{r')'P{r') 
h J 



‘^^V{r)d'{r) . 



(3.14) 



For the boundary conditions, we have 
|r ■ 



’ — r'|= \J — 2r ■ ‘ 



J I «/2 



= r\\- 



1 - 



2r • r' 



o\-. 



Note that in the last equation, 



r = — . 

r 



(3.15) 



(3.16) 



Let us introduce /C, the wave number in the direction of the outgoing unit 
vector f , 

K = Kr . (3.17) 

So the Green function for the operator expressed by (3.9), has the 

following asymptotic behaviour: 

exp{iKr — iK ■ r') 



g{r,r') 



47rr 



(3.18) 



Let us now introduce the asymptotic behaviour of the Green function 
(3.18) into the integral equation (3.8): 



d'{r) 



exp(ii<-2)-^ 



.( 3 . 19 ) 



47rr 



From the equation 

3 ,exp(iKr — iK, ■ r') 



d'^r' 






V{r')'^{r') 



r 47 t 



(3.20) 
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we can conclude that, if the scattering amplitude is given by 

J d^rexp{-iK.-r)V{r)<I'{r) , (3.21) 

then the boundary conditions are satisfied. In (3.21), /C is the wave number 
of the scattered particles, and If'(r) is the scattering wave function. 

Let us suppose that the ratio between the electron kinetic energy and the 
atomic potential energy is high enough to render the scattering weak and 
'f'(r) not very different from the incident plane wave exp(iKz). This is the 
assumption which is the basis of the first Born approximation, i.e. 

'f'(r) = exp{iKz) = exp{iK ■ r) . (3.22) 

Utilising the first Born approximation, expressed by (3.22), the previous equa- 
tion (3.21) becomes 

= [ d^rexp{—iK-r)V{r)exp{iK-r). (3.23) 

2ttU J 

If we use Tiq to indicate the momentum lost by the incident electron. 



(3.24) 

(3.25) 

(3.26) 



hq = h{K-K) , 

for fast particles we can write that 

f{d, (/>) = - -^2 J exp(iq • r) V (r) . 

As we are interested in a central potential, then 
U(r) = V(r) , 
and as a result, 

m<p)= m 

m 

= 2 d(j) sin 9 d9 / dr exp{iqr cos 9) V{r()3.27) 

2Trh jQ jQ Jo 

We carry out the integrations over (j) and over 9 and obtain 
2m 



f{9) = — / sin{qr)V{r)r dr . 

Ti q Jo 



(3.28) 



We are interested in the calculation of differential elastic scattering cross- 
sections in the first Born approximation for a screened Coulomb potential, 
such as a Wentzel-like potential [1], 



V\f) = exp I — 

r 



(-D 



(3.29) 



The exponential factor here represents a rough approximation of the screening 
of the nucleus by the orbital electrons, while the a parameter is 
ao 



a = 



where oq = /me^ is the Bohr radius. 



(3.30) 



3.1 The First Born Approximation 
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Let us calculate the scattering amplitude: 

f{S) = f sin{qr) exp(--) dr . 

n Q JQ \ a/ 

As the equation 



^ / r \ 

sin(gr) exp(^ j dr 



{1/ay 



holds, we can conclude that 






2^4 



= I f(0]\^ = — 

dl2 [q^ + {l/ayy 



(3.31) 



(3.32) 



(3.33) 



On the other hand, \K\ = \K\ and q = K — K, and, as a consequence 
q^= {K - K.) ■ {K - K.) 

+ -2KK.COS9 

= 2K^(l-cos6») , (3.34) 



where 9 is the scattering angle. 

The electron kinetic energy is given by 



2m 



(3.35) 



so that the differential elastic scattering cross-section for the collision of an 
electron beam with a Wentzel-like atomic potential is given in the first Born 
approximation by 



da 1 

df2 (1 — COS0 -I- a)^ 



(3.36) 



In the last equation, a is the screening parameter, given by 
1 me'^TT^ 

2K^ ^ ■ 

The well-known classical Rutherford formula, 

da Z^e^ 1 
dO ^ 4A2 (l-cos6»)2 ’ 



(3.37) 



(3.38) 



can be obtained by imposing a = 0 in (3.36). 

The total elastic scattering cross-section can be obtained from 



^el — 



/ 



>■ 



(3.39) 



For a Wentzel-like potential, the total elastic scattering cross-section can be 
easily calculated thus: 
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^ el 



AE^ 




sin § M 



1 

(1 — cost9 + ay 



TrZ'^e'^ 1 
£;2 a{2 + a) ' 



(3.40) 



If a ^ 0 and, as a consequence, the differential elastic scattering cross- 
section is given by the classical Rutherford formula, the total elastic scattering 
cross-section diverges, reflecting the long-range nature of the pure Coulomb 
potential. 

The elastic mean free path is the reciprocal of the total elastic scattering 
cross-section divided by the number JV of atoms per unit of volume in the 
target: 



1 _a(2 + a)E^ 

Naei Nne'^Z'^ 



(3.41) 



3.2 The Density Matrix and Spin Polarisation 

Following the previous section concerning the scattering of a spinless particle, 
it is necessary to introduce the density matrix, the spin polarisation and 
relativistic phenomena in order to take into account the effect of spin on the 
elastic scattering. 

The spin orientation of an electron beam is typically known through a 
probability distribution, which means it cannot be specified by a single state 
vector. Instead, it is a quantum system in a mixed state, constituted by n 
subsystems, each of which is in a pure state. 

Let us indicate by |a) the normalised state vectors of the systems in pure 
states. If we denote by |i) a complete set of orthonormal eigenvectors, then 
we can expand each pure state as follows: 

i 

where the coefficients of the expansion Ci are the projections of |a) on |i): 

Ci = (i|a) . (3.43) 

Let us now consider an operator A and carry out an evaluation of its average 
value 

n 

{A) = '^pa{a\A\a) , (3.44) 

a—1 

in which we have denoted by Pa the probability of obtaining the pure state 




3.2 The Density Matrix and Spin Polarisation 
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Let us now define the density operator 

n 

a—1 

Before proceeding, we can observe that the matrix elements of the density 
operator constitute the density matrix and are given by 

n n 

= y^(t|a)Pa(g|t') = '^PaC*Ci . (3.46) 

a=l a=l 

It is possible to see that the average value of operator A can be calculated 
as the trace of operator pA, once the density operator is introduced. The set 
of orthonormal eigenvectors denoted by |t) is supposed to be complete. As a 
consequence, 

^|*)(*| = 1. (3.47) 

i 

So, the average value of operator A is given by 

n 

(A) = '^Y^Y^{i\a)pa{a\i'){i'\A\i) 

a—1 i i' 
i i' 
i 

= Tr(^A) . (3.48) 

If we now consider the spin space, a complete set of 2 x 2 Hermitian 
matrices is given by the Pauli matrices Ci and the unit matrix I, so that the 
density matrix may be written as 

p = vqI + '^ Viai , (3.49) 

I 

where the Vi, with i = 0, 1, 2, 3, are real coefficients. 

The average values of the Pauli matrices are given by 

(CTi) = Tr(^CTi) . (3.50) 

The value of vq can be obtained by observing that 

Tr(/x) = Tr(Ml) = (/) = !, (3.51) 

and 

Tr(cTi) = 0 . (3.52) 

As a consequence, 

1 = Tr(^) = Wo Tr(I) = 2vq , 



(3.53) 
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or 

Vq=^- (3.54) 

Concerning the calculation of the values of the parameters Vi with i = 
1,2,3, we can observe that 

Tr(cr*cTj) = 2% . (3.55) 

Consequently 

(CTi) = Tr(^CTi) = 2vi , (3.56) 



or 



1 



Vi = ■ 



(3.57) 

If we now introduce the components of the polarisation vector, defined as 
P, = (a,) = Tr(/xa,) , (3.58) 

the density matrix may be written as 






or 



= 



1 ( I + P3 
2 I Pi + iP2 



Pi - lP2 
1 -P 3 



(3.59) 



(3.60) 



Choosing the z axis in the direction of the polarisation vector P, we obtain 
Pi = P 2 = 0, P 3 = |P| = P and 



^=2 



1 fl + P 
0 



0 

1-P 



The diagonal elements of the density matrix, 

Uti = , 



(3.61) 



(3.62) 



represent the probability that an electron of the beam is found in state i. As 
a consequence, if u indicates the number of electrons of the beam with spin 
up, and d the number of electrons of the beam with spin down, we can write 

1 + P u 



2 u + d ' 

so that the polarisation P of the beam is given by 
P=^~^ 



u + d 

Note that /r can be decomposed as follows: 

' 1/2 



(3.63) 



(3.64) 



Ai=(l-P)Mi/2 + /^m = (l-^^)(^^^ 0 ) • (3-65) 




3.3 Relativistic Partial- Wave Expansion 



25 



The diagonal elements of the matrix /r.1/2, 

correspond to a completely unpolarised system, since the probabilities of spin 
up and spin down are in this case equal to 1/2. The diagonal elements of the 
second matrix /ii, 

= (J 0) ' 

correspond to a pure state in which all the electrons of the beam have spin 
up; i.e. for a density matrix equal to /ii, the spin projection is in the direction 
of P for all the electrons of the beam. The beam in this case is said to be 
totally polarised. When P = 0 and the density matrix is ^1/2, the beam 
is completely unpolarised. When P = 1 and the density matrix is or 
P = —1, the beam is totally polarised. In all other situations, corresponding 
to 0 < |F| < 1, the beam is said to be partially polarised. 



3.3 Relativistic Partial- Wave Expansion 



The fundamental equation of relativistic quantum mechanics is the Dirac 
equation. The wave function, as is well known, is a four-component spinor, 
and the asymptotic forms of the four components of the scattered wave are 

~ aiexp{iKz) + bi{0,(l)) — (3.68) 

r 

r — >-oo 



The differential elastic scattering cross-section is given by 

da _ l&iP + I&2P + 
df2 |aip + |a2p + |o3p -I- |o4p 

_ \h\^ + I& 2 P + c|6ip -k c|&2p 

|aip + |fl2p + c|aip -k c|a2p 
I&1P+I&2P 
|aiP + |a2P ’ 



(3.69) 



where c is a constant of proportionality which takes into account the fact that 
the Oi and the bi coefficients are not all independent. Indeed, asymptotically 
the scattered wave is made up of plane waves proceeding, from the centre, in 
various directions; and the coefficients of the solutions for a plane wave are 
not all independent. 

If the spin is parallel to the direction of incidence (spin up), oi = 1, 
02 = 0, 61 = /+(0, (/>), 62 = g'^{9,(j)), where /+ and are two scattering 
amplitudes. 
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The asymptotic behaviour is described by the following equations: 

I'l ~ exp{iKz) + /+( 6 >, ^ (3 70 ) 

r — >oo 

r 

r — >-00 



The case of spin antiparallel to the direction of incidence (spin down) corre- 
sponds to fli = 0, 02 = 1, bi = g~{9,(j)), 62 = f~{9,cj)). The asymptotic 
behaviour is now given by 



ifi 



9 {9 A) 



exp{iKr) 

r 



(3.72) 






exp{iKz) + f {9,(j)) 



exp{iKr) 



(3.73) 



The Dirac equations for an electron in a central field are given by the following 
(see Chap. 2): 

1 -I- k 

[E + m-V {r)]Fi^ (r) -I ^ -I — Gf (r) = 0 , (3.74) 

-[E-m- V{r)]Gf{r) + \ ^ ’ + F±(r) = 0 . (3.75) 

dr r 



The natural units fi = c = 1 are used here, as they are particularly convenient 
for the quantum-relativistic equations. The superscript “-I-” refers to the 
electrons with spin up (fc = — (/ -I- 1)) while ” refers to electrons with spin 
down (k = 1). Once the new variables 



fj.{r) = E + m — V{r) , 
i/(r) = E — m — V(r) , 

, dn 

have been introduced, the Dirac equations become 



and 



\ dr r 



dEj^ _ n' f dGf ^ 1 -f 



dr 



dr 



+ 



1 f d^Gf 1 + kdGf l + k : 
V r dr j -2 i 



(3.76) 

(3.77) 

(3.78) 



(3.79) 



r dr 



(3.80) 
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Therefore, after simple algebraic manipulations, we obtain the following: 

dGf , f _ k{k + l) l + k^l'' 






+ I - - — 

V r /r 



dr 



liiy - 



A* 



^ Gf =0 .(3.81) 



Let us now introduce the function where 











1/2^1 



Upon observing that 
= E^-m^ , 
it is possible to see that 

Hi2 = - 2EV + . 

We conclude that, once the function U^{r) has been defined, i.e. 

-Ut(r)=-2EV + V--^-t^ + Y--Y4. 

‘ r 2 fjL i fj? 

the following equation holds: 

^ fc(fe -I- 1) 



dr 2 



+ K^-Uf{r) 



Gt = ^- 



(3.82) 

(3.83) 

(3.84) 

(3.85) 

(3.86) 



For large values of r, Qf^ is essentially sinusoidal. Indeed, when r is large 
enough, V{r) is negligible, is almost constant and the solution of the 
equation is therefore a linear combination of the regular and irregular spher- 
ical Bessel functions multiplied by Kr (see Appendix C). Taking account of 
the fact that G)*", we can therefore conclude that 

~ ji{Kr) cos — ni{Kr) sin -qf . (3.87) 



Here ? 7 ; are constants to be determined. Taking into account the asymptotic 
behaviour of the Bessel functions. 



ji{Kr) - ^ sin fiFr-^ 

r — KXD ^ 



ni{Kr) 



1 f 

cos I Kr 

Kr V 2 



(3.88) 

(3.89) 



we can conclude that 

Gf ~ sin (ATr — cos cos ( A"r — ^ ) sin ryf .(3.90) 

J\. T' \ ^ ^ \ " 
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Therefore, 



G 



+ 

i 



r — *oo 



sin I Kr 

Kr \ 




and 




r — KXD 



1 

— sin 
Kr 





(3.91) 



(3.92) 



The phase shifts represent the effect of the potential V (r) on the phases 
of the scattered waves. 

Before proceeding, we need to demonstrate the following equation: 

OO 

exp(iKrcosd) = ^^(2l + l)i‘ji(Kr)Pi(cos0) , (3.93) 

1=0 

where Pi{cos6) are the Legendre polynomials (see Appendix C) and ji{Kr) 
the spherical Bessel functions. In order to demonstrate this, let us first observe 
that a plane wave describing a free particle with the z axis in the direction 
of K may be expressed as an expansion in a series of Legendre polynomials 
Pi{cos9): 

OO 

exp(*ATz) = exp{iKr cos 0) = ''^ciji{Kr)Pi{cosO) . (3.94) 

1=0 

Let us define the two new variables s = Kr and t = cos 9, to have 

exp{ist) = '^ciji{s)Pi{t) . (3.95) 

i 

Differentiating with respect to s, 

itexp{ist) = itciji{a)Pi{t) = c; P;(t) . (3.96) 



Recalling one of the properties of the Legendre polynomials (see Appendix 
C), we have: 



m) = 

Therefore, 



{I + l)Pi^i(t) + IPi-iit) 



t{2l + l) 



(3.97) 



itexp{ist) = E itciji{s) 

i 



(I + l)Pi+i(t) + lPi-i(t) 

t(2/ + l) 

I ■ / \ /+1 

Cl-ijl-i(s) + — — Q+iJZ+i(s) 



2/-1 



21 + 3 



. (3.98) 
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On the other hand, it is well known that 

dji{s) I . 1+1 . , . 

- - — — ji+i(s) , 



ds 21+1 
and as a consequence, 



itexp{ist) = ciPi{t) 



21+1 

I 

21 + 1 



! \ ^ + 1 ■ / ^ 
J;-i(s) - ^^x + l^i++^> 



(3.99) 



(3.100) 



Consequently, from (3.98) and (3.100), the following is obtained: 



Haw 






Cl 



+ 1 ) 



lCl-1 

21 + 1 21-1 

Cl ici+i 



21 + 1 2/ -k 3 



= 0 . 



(3.101) 



The Legendre polynomials Pi (t) are linearly independent (orthonormal) , and 
therefore 






Cl 



— jl+l(s) {I + 1) 



^Cl-l 

21 + 1 21-1 
Cl 



*Q+i 



^2/-l-l 2l + ?>, 

Every value of s satisfies the last equation if and only if 



1 



-Cl = 



-Cl-l . 



(3.102) 



(3.103) 



21+1 21-1 

In order to obtain an explicit expression for c; we need to know the value of 
the first coefficient of the set, i.e. cq. Imposing r = 0 in (3.94), we obtain 



exp(O) = 1 = ^ciji(Q)Pi{cos,e) . 



(3.104) 



Since ji{Q) = 0 for any / yf 0, while jo(0) = 1 and Pq(cos0) = 1, we may 
conclude that cq = 1. The recursive repetition of the relation (3.103) allows 
us to obtain the values of the coefficients. 



Cl = { 2 l + iy , 

and the expansion of the plane wave in Legendre polynomials. 



(3.105) 



exp(*prrcos0) = exp{iKz) = ^(2Z + l)i''ji{Kr) Pi {cos 9) . (3.106) 



1=0 



Let us remind the reader that we are looking for the functions <P'i and 
d '2 which satisfy the asymptotic conditions. So, we must begin by expanding 
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them in spherical harmonics (see Appendix C): 

C30 

+ BiGT]Pi{cose) , 



Z=0 

oo 



^2 = ^[QG+ + AGr]P;'(cos0) exp(i^) . 



(3.107) 



(3.108) 



The coefficients Ai, Bi, Ci and Di can be determined by considering the 
asymptotic behaviours of the functions involved. Let us begin with the func- 
tion and observe that 



<I'i - exp(iKz) = '^[AiGl + BiGi - {21 + l)i‘ ji{Kr)] Pi {cos 0) . 

(3.109) 



1=0 



As 



'I'l — exp{iKz) 



exp{iKr) ^ 



r{0.^) 



(3.110) 



the following occurs: 



1 

Kr 



E 

;=o 



hr 



Ai sin [ Kr - — + rj'f \ + Bi sin [ Kr - — + rji 



Itt 



— {21 + l)i^ sin ( Kr — — 



Pi{cos9) 



exp{iKr) ^ 



r(^,</>) • 



As a consequence, 

exp(f ATr) 
2iKr 



OO j 

^exp(-iy) 

i=0 



X [Ai exp{ir]^) + Bi exp{ir]i ) - {21 + l)i']P;(cos 6) 
exp{—iKr) hr. 

1=0 

X [Ai exp{-ir]^) + Bi exp{-irj{^) - {21 + 1) i’'] Pi {cos 9) 
-/+(0,<(>) • 



exp(iA>) = 



The asymptotic conditions are satisfied if 

Ai exp(-i? 7 ;+) -I- Bi exp{-ir]'f) = {21 + l)i‘ . 



(3.111) 



(3.112) 



(3.113) 
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With the choices 

Ai = li’' exp{irj'^) , (3.114) 

B[ = {I + l)i'‘ exp{irif ) , (3.115) 

(3.113) is satisfied. 

Proceeding in a similar way regarding the W 2 function, we can therefore 
choose 

C; = i'exp(z? 7 ;+) , (3.116) 

Di = —i^exp{irif) . (3.117) 

In conclusion, for electrons with spins parallel to the direction of incidence, 
we have 

00 

'f'l = + 1) exp(i77;“)G)" -I- lexp{iri'^)G'^]i’‘Pi{cosO) , (3.118) 

1=0 

00 

'f '2 = y~^[exp(t? 7 )^)G)^ — exp(z?7;“)G)"]i*P/(cos0) exp(i</>) , (3.119) 



1=1 

and, by using (3.112), 

/+(0,<^) =f+{6) 

1 °° 

= ^ XI + l)[exp( 2 i? 7 ") - 1] 

+l[exp(2ir]'^) — l]}P;(cos0) , (3.120) 

. C50 

9^ {&,<!}) = -exp{2ir]^)]Pl{cos0)exp{i<j)) . (3.121) 

1=1 

For electrons with spins antiparallel to the direction of incidence (spin down), 
where we indicate the scattering amplitudes by f~ and g~ , it is possible to 
see that 

r(0,0) = /+(0,<^) (3.122) 

and 

9~{0,(p) = - 9 ~'~ (0,(1)) exp . (3.123) 

It is therefore convenient to define the functions 

00 

f{e) = Y,M{cose) , (3.124) 

1^0 
00 

5(0) = ^P;Pi'(cos0) , 

1=0 



(3.125) 
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where 



+ l)[exp( 2 t 77 ; ) - 1 ] + /[exp( 2 i? 7 +) - 1 ]} 

I3i = ^{exp(2ir?+) - exp(2i77")} . 

With this notation, we have 

t = r = f. 

=g exp(t 0 ) 

and 

g~ = ■ 



(3.126) 

(3.127) 

(3.128) 

(3.129) 

(3.130) 



For an arbitrary spin direction, the electron incident plane wave will be 
given by 'Pi = Aexp{iKz) and P 2 = Bexp{iKz), and as a consequence 
oi = A, a 2 = B. Furthermore, 



hi = Af+ + Bg = Af - Bg exp(-i(/>) 
&2 = Ag+ + Bf~ = Bf + Agexp{i(l)) . 

Consequently, 

da 

dfl 

= {\f\^ + \g?)[l + iS{9) 



(3.131) 

(3.132) 



AB* exp{i(j)) — A* B exp(— i(/>) 






\B\ 



where S{9) is the Sherman function, defined by 

j9*-f*g 



S{9) = i 






Note that 



. AB* exp(i(/>) — A* B exp{—icj)) 



= ^\a 2 cos 4> — ai sin (ji)^ 



3 . 133 ) 



(3.134) 



(3.135) 



where cti, a 2 and (T 3 are the Pauli matrices and ^ is the two-component spinor 

(3.136) 






B' 



e = 



B/V\W- 

A* 

t/\^TW7\WW\ 



B* 



(3.137) 



As the 2 ; axis has been chosen along the incidence direction, the unit vector 
perpendicular to the plane of scattering is given by 

n = (— sin (j), cos <j), 0 ) , 



(3.138) 
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so we can therefore write 



((72 cos 4> — ai sin = P fi 



( 3 . 139 ) 



where P is the initial polarisation vector of the electron beam. The differential 
elastic scattering cross-section can then be recast in the following form: 

^ = i\f\^ + \gni + S{e)P-n]. ( 3 . 140 ) 

Note that, if the beam is completely unpolarised, then P = 0 and 

^ = I/I^ + I.P. ( 3 . 141 ) 

The total elastic scattering cross-section (aei) and the transport cross-section 
(atr) are defined by 



(Tel — 27T 



■ sin 6 d6 , 



/*7T 

atr = 2 tt / (1 — COS 0 )— — sin 0 d 0 , 

Jo dS 2 



( 3 . 142 ) 



( 3 . 143 ) 



which can be easily calculated by numerical integration. 

Note that by imposing 

dT=vt = Vi ( 3 . 144 ) 

in the previous equations, we can obtain the non-relativistic results. Indeed, 



+ l)[exp(2ir?i) - 1] + l[exp{2i-qi) - 1]} 

2iK 



= ^(2^ + l)[exp(2i7?i) - 1] 



Bi = 0 



( 3 . 145 ) 

( 3 . 146 ) 



so that 



f{d) = ^ + l)[exp{2ir]i) - l]Pi{cos0) 

* 1=0 
OO 

= —'^{21 + 1) exp{ir]i) sim]iPi{cos0) , 



( 3 . 147 ) 



g{0) = 0 



( 3 . 148 ) 



( 3 . 149 ) 
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3.4 Calculation of the Phase Shifts 



In order to proceed, it is convenient to perform the following transformation 

[ 2 ]: 



F±(r) = «f(r) 



sin (r) 



Gf(r)=af(r) 



cos (r) 



After simple algebraic manipulations, (3.74) and (3.75) become 

1 daf(r) 



[E + m— V(r)] tan (pf (r) 



af{r) dr 



_tan</)f(r)^^^ 



^ = 0 
r 



(3.150) 

(3.151) 



(3.152) 



— [E— m — V{r)] cot 4>^{r) + 



+ cot 



dr 




1 daf{r) 
af{r) dr 



and therefore 
dr 



k 

r 



sin2^)*^(r) 



m cos 2(j)f- (r) + E — V (r) , 



1 daf{r) 
af{r) dr 



k 

r 



cos2^)*^(r) 



msin2(j}^{r) . 



(3.153) 



(3.154) 

(3.155) 



For 0 < r < fi/mc, the spherical symmetric electrostatic potential experi- 
enced by a point charge at distance r from the nucleus, P(r), may be ap- 
proximated by the following: 



V{r) ~ 

>^0 



Zq + Zir + Z-ir"^ + Zsr^ 
r 



(3.156) 



Expressing the electrostatic potential as the product of the potential of a bare 
nucleus multiplied by a screening function ^(r) having the analytical form 

p 

= exp(-«ir) , (3.157) 



^ Ai — 1 , 
i=l 



(3.158) 
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we can easily evaluate Zq, Zi, Z^ and Z 3 : 




Zq = Ze^ ^^Ai = Z(? , 

i 


(3.159) 


P 

Zi = —Zq aiAi , 


(3.160) 


^2 = , 

2=1 


(3.161) 


2 » = -f 

2 = 1 


(3.162) 


Let us expand (j)f^ as a power series 




(’’) = </>w + + . . . . 


(3.163) 


It is possible to see, after simple algebraic manipulations 
ships between the coefficients of this expansion and Zq, 
the following [3]: 


, that the relation- 
al, Z 2 and Zq are 


sin 2 (/),o = - y , 


(3.164) 


±. E + Zi — m cos 

1 — 2kcos2(j)^ ’ 


(3.165) 


± 2(j)f^ sin 2 </^(m - + Z 2 

2 — 2kcos2(jy^ 


(3.166) 



A± 20^sin2(/)^(77i — 2fc(/)^) + 20^^ cos 2</)^[m — (2/3) 

^13 ~ 



3 — 2fccos2(/))^ 



with the extra conditions 

0 < 2</± < ^TT 

if fc < 0, and 

7T < 2(j)% < ^7T 

if fc > 0. 

Let us now calculate the phase shifts, examining (3.151): 



g;± = 



cos(/f(r) af ± ,± af cos(/f(r) 



- 4-sin</); (r)</; (r) - 



(3.167) 

(3.168) 

(3.169) 

(3.170) 
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SO that 

^ {r)~ - , (3.171) 

or 

= —{E + m — V) tan (j)f{r) — — . (3.172) 

G, r 

Let us observe that the asymptotic form of the solution in the regions corre- 
sponding to large values of r for which F(r-) « 0 is (see (3.87)) 



Gf" = ji{Kr) cos-qf" — ni{Kr) sinqf- , (3.173) 

where — m^, are the /th phase shifts, and ji and ni are respec- 

tively the regular and irregular spherical Bessel functions. Therefore, 



G)^ Kj'i (Kr) cos qf" — Kn[ {Kr) sin qf 

Gf ji{Kr) cos qf" — ni{Kr) sin qf- 


(3.174) 


Taking into account the properties of the Bessel functions 




3'i{x) = -ji{x) - ji+i{x) , 

X 


(3.175) 


n'i{x) = -ni{x) - ni+i{x) , 

X 


(3.176) 


we may conclude that 




tann± - Kji+,{Kr) - UKr){G[^ /G^) 

(l/r)ni{Kr) - Km+,{Kr) - ni{Kr){G[^ /Gf) ' 


(3.177) 


Let us define 




= lim . 

r — >oo 


(3.178) 


For large values of r, (3.172) becomes 




G(^ / X 7-1- 1 + 

± = {E+m)tan<j)^ 

(_r, r 


(3.179) 



i 



and therefore 

± Kji+i{Kr) - ji{Kr)[{E + m) tan^f + {1 + I + k)/r] 

tan 71, = • 

Kni+i{Kr) - m{Kr)[{E + m) tan(j)f -|- (1 -|- Z -|- k)/r] 

(3.180) 

Using this last equation, we can calculate the phase shifts of the scattered 
wave and, therefore, the functions f{9), g{9) and the differential elastic scat- 
tering cross-section. 
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3.5 Exchange and Solid State Effects 

The incident electron may be absorbed by the atom while another electron 
is emitted, so that such an exchange effect modifies the wave equations. This 
effect may be approximated, as proposed by Furness and McCarthy [4], by 
adding the following to the electron-atom potential energy: 

K. = - Vs) -^[{E- Vs)^ + 47 Tpe %]^/2 ^ 

Here E represents the electron energy. Vs the electrostatic scalar potential, p 
the atomic electron density, e the electron charge and oq the Bohr radius. 

Solid state effects should also be taken into account when the target atom 
is bound in a solid. An approximate way to treat this effect is to introduce 
what is known as central muffin-tin model, in which the potential of each 
atom of the solid is altered by the nearest neighbours; assuming that they 
are located at a distance 2rws, the resulting potential is 

VsoHdir) = V(r) + V(2rws ~ r) - 2V(rws) , (3.182) 

for r < rws, and zero elsewhere. Here rws is the Wigner-Seitz radius, 

rws = 0.7346(A/5)^/3^ (3.183) 

where A is the atomic weight and <5 the mass density expressed in g/cm^. 
The term 2V(rivs) also has to be subtracted from the kinetic energy of the 
incident electron. 



3.6 Comparing Theory and Experimental Data 

In Figs. 3. 1-3.5 we show a comparison of the electron differential elastic 
scattering cross-section calculated both from the screened Rutherford formula 
(first Born approximation) and by the relativistic partial-wave expansion 
method, with the experimental data available. The data concern different 
electron energies, target atoms and molecules. The energies considered are 
all below the Born threshold and it is evident from the comparison with 
experiment that the relativistic partial-wave expansion method allows a much 
more accurate description of the elastic scattering for these energies than does 
the screened Rutherford formula. 

On the basis of this and other comparisons, and taking into account the 
experimental errors, we can say that the accuracy of the differential elastic 
scattering cross-sections is expected to be 1-2%, for scattering angles greater 
than 5°, as calculated by the relativistic partial- wave expansion method. 
The effect of the charge-induced atomic polarisation was neglected in our 
calculation code: such an effect may be relevant for angles lower than ~5° 
and this means that we can expect an accuracy in the total cross-section of the 
order of 5-6%. In Fig. 3.6, we show a comparison of the total elastic scattering 
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Fig. 3.1. Differential elastic scattering cross-section dajdQ of 400 eV electrons 
scattered by Ar as a function of the scattering angle Q. Dotted line: screened Ruther- 
ford formula. Solid line: relativistic partial- wave expansion method [5]. •: Iga et al. 
experimental data [6] 




Fig. 3.2. Differential elastic scattering cross-section dajdfi of 1000 eV electrons 
scattered by Ar as a function of the scattering angle 9. Dotted line: screened Ruther- 
ford formula. Solid line: relativistic partial- wave expansion method [5]. •: Iga et al. 
experimental data [6] 
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Fig. 3.3. Differential elastic scattering cross-section dajdfi of 1100 eV electrons 
scattered by An as a function of the scattering angle 9. Dotted line: screened Ruther- 
ford formula. Solid line: relativistic partial- wave expansion method [.5]. •: Reichert 
experimeirtal data [7] 




Fig. 3.4. Differential elastic scattering cross-section dajdfi of 300 eV electrons 
scattered by Hg as a function of the scattering angle 9. Dotted line: screened Ruther- 
ford formula. Solid line: relativistic partial- wave expansion method [5]. •: Bromberg 
experimental data [8]. o: Holtkamp et al. experimental data [9] 
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Fig. 3.5. Differential elastic scattering cross-section da/dfi of 500 eV electrons 
scattered by CO as a function of the scattering angle 9. Dotted line-, screened 
Rutherford formula. Solid line-, relativistic partial-wave expansion method [5]. •: 
Bromberg experimental data [10] 




Fig. 3.6. Total elastic scattering cross-section a^i of electrons scattered by Ar 
as a function of the electron kinetic energy E. Solid line-, relativistic partial-wave 
expansion method [5]. -|-: Jansen et al. experimental data [11], •: DuBois and Rudd 
experimental data [12]. o; Iga et al. experimental data [6] 
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d (deg) 

Fig. 3.7. Sherman function S(Q') for 1500 eV electrons scattered by gold as a 
function of the scattering angle Q. Solid line: relativistic partial-wave expansion 
method [5]. •: Deichsel and Reichert experimental data reported by Bunyan and 
Schonfelder [3] 



cross-section calculated using the relativistic partial-wave expansion method 
and the experimental data of Jansen et al. [11], DuBois and Rudd [12] and 
Iga et al. [6]. 

The first transport cross-section is very important for multiple-scattering 
processes: this transport cross-section does not appear to be very much in- 
fluenced by the inaccuracies of the differential elastic scattering cross-section 
at low scattering angle (owing to the factor (1 — cos0) in (3.143)), which are 
a consequence of neglecting atomic polarisation. 

We are able to take into account the polarisation phenomena with the 
relativistic partial-wave expansion method. In Fig. 3.7, we compare the Sher- 
man function (calculated with this method) with the Deichsel and Reichert 
experimental data reported by Bunyan and Schonfelder [3] . 
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4 Inelastic Scattering 



4.1 The Classical Theory 



Assume that a free electron is at rest. At a distance r from it, another electron 
is travelling along the 0 direction. If 9 represents the angle between r and z, 
V represents the velocity of the incident electron and the instant of impact is 
t = 0, 



— vt = r cos 9 . 

Defining the impact parameter b by 
h = r sin 9 , 



(4.1) 

(4.2) 



it is possible to see that b is the distance from the trajectory of the incident 
electron to the target: 

7,2 I 2j.2 2 

0 + V t = r . 



F indicates the force of repulsion between the two electrons, where 



F. = 



‘ cos 9 



and 



F^, = 



sin 9 



(4.3) 



(4.4) 



(4.5) 



and Pz and Px are the components of the momentum transferred to the target 
electron. As 



_ vt 

cot 0 = ^ 

b 


(4.6) 


differentiating with respect to 9, we can obtain 




dt b 


(4.7) 


d9 V sin^ 9 


Therefore, 




1 dt 1 


(4.8) 


d9 bv 
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The X component of p is thus given by 

„2 2e^ 



Px = 



Fr dt = 



bv , 



sin OdO = 



0 



bv ’ 



(4.9) 



while the 2 : component of p is zero. The energy W transferred to the electron 
at rest is 



W = 



(4.10) 



Px e 
2m b'^E ’ 

where E = mv^ j 2 is the kinetic energy of the incident electron. 

The number of electrons in the volume 2'Kb db dz is given by 2i:b db dz N Z, 
where N is the number of atoms per unit of volume in the target and Z the 
target atomic number. Thus, the energy lost by the incident electron per unit 
of length in the target (called the stopping power) may be calculated from 

2Tre'^NZ 



rl C 

— = / 2TibdbWNZ = 

dz I 



E 



T 



(4.11) 



where &min and 6max are, respectively, the minimum and the maximum impact 
parameter. 

The atomic electrons may be considered as a set of oscillators with fre- 
quencies Vi and amplitudes fi. As a consequence, we can rewrite the last 
equation as 



dE 2'ire'^N 



dz 



E 



E/* 



~b 



(4.12) 



The sum of the amplitudes of the oscillators has to be equal to the atomic 
number Z\ 



E/* = ^ 



If hvi indicates the binding energy of the tth atomic electron, then 

^max = ■ 

The minimum impact parameter is given by 
&min = h/mv , 
so that 



dE 2'ire'^N 



E 



mv 

hvi 



dz E 

Let us define the mean ionisation potential as 



(4.13) 

(4.14) 

(4.15) 

(4.16) 



(4.17) 




4.2 Dielectric Function and Stopping Power 



45 



SO that the stopping power becomes 



dz E 



dE 2-Ke^NZ , 
in 




(4.18) 



This equation, deduced here with the classical theory, is substantially 
correct. Using a quantum mechanical treatment, the equation becomes the 
following: 



where the constant C equals «1.166 (Bethe equation) [1]. The Bethe formula 
is valid for energies higher than ~ J, as it reaches a maximum for E « 2.3J, 
and then goes to zero for E k, J/ 1.1 66: for lower energies the predicted 
stopping power becomes negative. Therefore, the low-energy stopping power 
requires a different approach, and its evaluation is based on the optical data 
of the material of interest. 



4.2 Dielectric Function and Stopping Power 

The response of a medium to an energy transfer fiuj and momentum transfer 
q is described by the complex dielectric function e{q,1luj), which is related 
to the optical data. If p{E, T%u>) is the probability for an energy loss tiu) per 
unit distance travelled by an electron of energy E, the stopping power is then 
given by the following equation: 



where the integration is over the allowed values of the energy transfer fiuj. If 
Tiuj is the energy transfer, q the momentum transfer and e(q, fko) the complex 
dieletric function, then p{E^ Tioj) is given by 




(4.19) 




(4.20) 




(4.21) 



where 




(4.22) 



and 




(4.23) 



Let us observe that 

e(0, flUj) = £i -I- 1£2 



(4.24) 
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where 



£i = Re[£r(0, huj)] = 

£2 = lm[£(0, huj)] = 2nk , 



(4.25) 

(4.26) 



and where n is the coefficient of refraction and k the coefficient of extinction. 
Therefore, 



Im 



-1 



e(0, hoj) 



lm[e(0, hiu)] 



£2 



2nk 



(4.27) 



|e(0,?icj)p ef + e| {v? + k'^Y 

Ashley [2, 3] showed that the relation between the dielectric function and 
the optical data may be approximated by 



Im 



-1 



_e{q, huj)_ 

pOO 

/ d{Tiuj') huj' Im 

Jo 



-1 


A 




_£r(0, huj')_ 




\ 2m /_ 



fiuj 



(.4.28) 



As a consequence, for E > 4Ef (where Ef is the Fermi energy), the stopping 
power may be calculated from 
^.2 rE/2 



dE 

ds 



■kTi E Jo 



Im 



-1 



Gr ( ^ ) h,u>d{hu>) . 



_e(0, fiuj) _ 

Ashley also gave an accurate approximation for the function Ge{x): 
Ge(x) 



(4.29) 



= In 



1.166 



X 



3 X /4 
- -X - - In - 

4 4 V a: 



ix3G-^lnflV-^-30) 

2 16 U/ 48^ ’ 



4.3 Inelastic Mean Free Path 



The electron inelastic mean free path can be calculated from the following 
equation: 

1 



^inel — 



(4.31) 



/ p{E, Tiuj) d{Tiuj) ’ 

where the integration is extended over all the allowed values of the energy 
transfer T%lo. In a way similar to that used to calculate the stopping power, it 
is possible to show that 



_ 

^inel 



2nh‘‘E Jo 



Im 



-1 



ie ( ^ ) d{hw) 



_£(0, huj) 

where the approximate evaluation of Le(x) proposed by Ashley is 

4\ 7 Q/n 33 

-X + x -^'^ - —0 

, a: y 4 32 



Lf.{x) = (1 — x) In ( — ) — ^x + a;3G — —x^ 



(4.32) 



(4.33) 
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Note that, for the calculation of the inelastic cross-section of low-energy elec- 
trons, Powell [4] proposes the following differential cross-section: 



dX, 



inel 



d{Tiuj) 



2ttTi^E 



Im 



-1 



e(0, huj) 



In 



cE\ 

Uto J 



(4.34) 



where c is a parameter depending on the energy loss Tilu = AE. 



4.4 Positrons 



The stopping power and the inelastic mean free path for positrons may be 
calculated in a similar way [2, 3] (except for the exchange). For 4,Ej < E < 10 
keV, 




p 



me^ 

2TTTi^E 




-1 

£t(0, huj) 



Gp 




Tiujd(Jiuj) , 



(4.35) 



i^inel')p 



r-E/2 



2tt?i^E 



Im 



-1 



e(0, hoj) 









(4.36) 



where 



and 



Gp{x) = In 



Lp{x) = In 



1 — X + 'Jl — 2x 
1 — X — Vl — 2a: 



1 — x/2 + y/1 — 2x 
1 — x/2 — ^Jl — 2x 



(4.37) 



(4.38) 



4.5 Plasma Oscillations 

As plasmon excitation is very important in electron energy loss, we would 
like to have a closer look at this phenomenon. Let us give a brief treatment 
of the problem, considering the gas of the valence electrons in a metal: if n is 
the number of electrons per cubic centimetre, then the collective excitations 
of this gas produce oscillations as a result of the passage of electrons. 

If ^ is a fluctuation of the distance r of the free-electron gas from the 
lattice of the positive charges of the solid, then the number v of electrons in 
a spherical shell of “thickness” ^ is given by the following, 

V = 47rr^n^ , (4.39) 

so that the retarding force is 

g2 

F = — -u = —A-Ke^n^ . 



(4.40) 
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As a consequence, 



m—k + 47re^n^ = 0 
dt^ 



(4.41) 



and hence the motion is a harmonic oscillation of frequency ujp, given by 



47re^n 



(4.42) 



A quantum of plasma oscillation is a plasmon. The energy of a plasmon is 
Ep = fiujp. Plasma oscillations confined to a metal surface have a frequency 
Ws, given by 

(4.43) 



LOs = 



V2' 



The mean free path for plasmon emission can be calculated by using the 
Quinn formula [5], 

-1 



Xp — 



2uq£' 



In 



LUqj 



{pj + 2mwp)^/^ -pf 
p — {p'^ — 2mu!py/'^ 



where 

p = {2mEy^'^, 

Pf = 

and Ef is the Fermi energy. 



(4.44) 

(4.45) 

(4.46) 



4.6 Comparing Theory and Experimental Data 



The theoretical results that we present in Figs. 4. 1-4.5 were obtained by 
following the method proposed by Ashley [2, 3] and using the mass absorption 
coefficients reported by Henke et al. [6]. 

The calculations of the refraction index n and of the extinction coefficient 
k were performed by use of [6] 



n = 1 — 



2TTmc^ 



X^NY^Xphp, 



(4.47) 



k = 



2'Kmc^ 



X^NY^Xphp. 



(4.48) 



In these equations, c is the speed of light, N the number of molecules per 
unit volume (each having Xp atoms) and A the photon wavelength, fip and 
f 2 p are, respectively, the real and the imaginary component of the atomic 
scattering factor. The calculation of the real and imaginary component of 
the atomic scattering factor was performed using the following equations [6] : 
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E (eV) 

Fig. 4.1. Inelastic mean free path Ai„e( of electrons in Si02 as a function of the 
electron kinetic energy E. +: theoretical calculation based on the optical data [7]. 
• : Flitsch and Raider experimental data [8] 




E (eV) 

Fig. 4.2. Inelastic mean free path Xinei of electrons in AI2O3 as a function of the 
electron kinetic energy E. +: theoretical calculation based on the optical data [7]. 
• : Battye et al. experimental data [9] 
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E (eV) 

Fig. 4.3. Electron stopping power —dE/ds in MgO as a fnnction of the electron 
kinetic energy E. +: theoretical calculation based on the optical data [7]. •: exper- 
imental data reported by Akkerman et al. [10], taken from Joy’s database [11] 




E (eV) 

Fig. 4.4. Electron stopping power —dE/ds in Si02 as a function of the electron 
kinetic energy E. -I-: theoretical calculation based on the optical data [7]. •: exper- 
imental data reported by Akkerman et al. [10], taken from Joy’s database [11] 
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Fig. 4.5. Electron stopping power —dE/ds in AI 2 O 3 as a function of the electron 
kinetic energy E. +: theoretical calculation based on the optical data [7]. •: ex- 
perimental data reported by Akkerman et al. [10], taken from the Joy’s database 



In these equations Ep is the incident photon energy, Z the atomic number, 
A the atomic weight, Na the Avogadro number and /i the photoabsorption 
cross-section (cm^/g). 

The Henke et al. data are given for discrete photon energies (beginning 
at 30.5 eV, and then 49.3 eV, 72.4 eV and so on), and we thus performed the 
calculation of the optical loss function utilising a cubic spline interpolation of 
the mass absorption coefficients. We performed polynomial extrapolation of 
the data to calculate the mass absorption coefficients for energies lower than 
30.5 eV [7]. 

At high energy, the stopping power is well described by the Bethe formula. 
With I = 142 eV for Si 02 , the stopping power at 10 keV obtained with the 
Bethe formula is 0.458 eV/A. The stopping power at 10 keV computed in the 
way described (i.e. using the optical data) is 0.462 eV/A, 0.9% higher than 
Bethe’s value. 

The inaccuracy of the calculated inelastic mean free path and stopping 
power is less than ~ 5 — 15%. 



[ 11 ] 




(4.49) 




(4.50) 
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5 Electrons Impinging on Solid Targets 



5.1 Backscattered Electrons 

The backscattering coefhcient is the fraction of electrons of the primary beam 
that emerge from the surface of an electron-irradiated “bulk” target, where 
the meaning of the word “bulk” is explained in the following equations. The 
mean path length travelled by an electron in a solid target when its energy 
reduces from the primary energy Eq to E is given by 



This equation allows us to calculate the maximum range of penetration, 



where ifmin is the energy at which an electron can be considered absorbed 
by the target for any practical purpose. 

When the thickness of the target is greater than R{Eq), it is a bulk target 
for the primary electron energy considered. In general, when the thickness is 
lower than R(Eq), the target is a film and the primary beam impinging on it 
is divided into three fractions, i.e. the backscattered, the absorbed and the 
transmitted electrons. 

If the target is a bulk target, then the number of electrons transmitted 
is zero, while the fraction of electrons backscattered reaches its maximum 
value (the backscattering coefficient). The remaining primary particles are 
absorbed, so that, if r is the backscattering coefficient, the fraction of ab- 
sorbed electrons equals 1 — r. 

Various semi-empirical approaches were proposed in the scientific liter- 
ature in the early 1960s to treat the problem of electron backscattering. In 
1960, Everhart [1] showed that the backscattering coefficient, for normal inci- 
dence, for high electron energy (greater than ~10 keV) and for low-medium 
target atomic number (lower than ~40-45), could be calculated from the 
following: 
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(5.1) 




(5.2) 
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r = 



g- 1 + 0.5° 
g + 1 



(5.3) 



where a = 0.045Z and Z is the target atomic number. 

The Everhart deduction is based on the semi-empirical Thomson- Whid- 
dington energy loss relation [2] 

El-E^ = cz , (5.4) 



where E represents the most probable electron beam energy at depth 2; 
in the target and c is a constant depending on the target. The Thomson- 
Whiddington energy loss relation can be obtained by assuming that the mean 
ionisation energy J depends linearly on E. The Bethe formula in such a case 
becomes 



dE k 
dz E 



where fc is a constant and, as a consequence, 



EdE 



Eq 



3 . 

2k 



E^ 



(5.5) 



(5.6) 



With c = 2k, this is the Thomson- Whiddington formula. Assuming that 
ddmin — 65 then 



Eq = cR 



(5.7) 



and, as a consequence, 
E^ = c{R - z) . 



(5.8) 



Let us introduce the supplement of the deflection angle, d = tt — 6. Tak- 
ing into account the Thomson- Whiddington equation, the Rutherford cross- 
section for a pure Coulomb field, da/df2, given by 



da 1 


(5.9) 


dn 4E2 (l-cos 6»)2 ’ 


becomes 




1 27t sin di? 


(5.10) 


^ 4c i? — z (1 + cosd)2 


The fractional range is defined as 




z 


(5.11) 



and the incremental number of electrons at a fractional depth u deflected 
through d is given by 

dv(d,u) = dN — da . (5.12) 

Here, N is the number of atoms per unit of volume in the target, so that 
dN = NS dz gives the number of atoms in the volume S dz, while v/S gives 
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the number of electrons incident on the surface, of area S. Once we have 
defined the quantity 

nZ'^e^N 



a = 



4c 



(5.13) 



the incremental number of electrons at a fractional depth u deflected through 
'd becomes 



v{u)du svadM 

dv{-d,u) = 2a- — ^ . 

1 — u (1 + cosd)"^ 



(5.14) 



On the plane at the fractional depth u, the number i>{u) of incident electrons 
can easily be calculated. Indeed, 

v{u) = V[)— / / diy{-d,u') , (5.15) 

Jo Jo 



and therefore 



r 

v[u) = vq — a 

Jo 



v{u')du' 



1 - It' 

The solution of this integral equation is as follows: 
v{u) = vo{l - u)“ , 
and hence 

dvid.u) „ 1 . svaddd 

^ ’ ’ = 2a(l - u)“"^ du ■ 



(5.16) 

(5.17) 

(5.18) 



z/Q (l + cosz?)2 

The total distance travelled by an electron reflected at an angle d at depth 
z into the target is s = z + z sec d. If dg is the maximum escape angle for a 
given z, then z(l + sec do) = E, or u(l + sec do) = 1. With do satisfying the 
relationship 

m( 1 + sec do) = 1 , 

the backscattering coefficient r is given by the following: 



r = 



J Vo Jo Jo 



° sin d dd 
0 (l + cosd)2 



(5.19) 

(5.20) 



The integration over d gives 



“ sin d dd 
and, hence 



1 



1 _ 1 

0 (1 + cosd)^ 1+ cos do 2 2^ 



/•1/2 /I 

2 a I (l-„r>(- 



— u ] du = 



g-l + (0.5)° 
g T 1 



(5.21) 



(5.22) 



The simple Everhart theory can be used to obtain further useful infor- 
mation. In 1978, Jablonski [3] showed that it can be utilised to calculate the 
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energy and angular distribution of backscattered electrons. As an electron 
reflected at depth 2 ; at an angle -d travels along a path given by 

s = Ru{l + secd) (5.23) 



and, in the Thomson- Whiddington approximation, 
25 



s = R 



1 - 



(5.24) 



then 



cos-d 
1 + cos 




(5.25) 



As a consequence, 

^ 2 cos -d E dE 

l+cosr9“^- 
If we define the function 




(5.26) 



(5.27) 



the energy and angular distribution of the backscattered electrons is then 
given by the following: 



d?v{'d,E) AaEi/Q / cost? p ^ sin -d cos -d 

dE dd Eq ■ 1+cosr?^^ ’ °7 (I+COS 73 



(5.28) 



5.2 Electrons in thin films 

5.2.1 Definitions, Symbols, Properties 

Let us consider an electron beam impinging on a solid target. If R represents 
the maximum penetration range, the target is considered a film if its thickness 
s is smaller than R. If s is greater than R, the target is a bulk target. The 
fraction of particles transmitted in such a case is zero, while the fraction 
of particles backscattered is equal to the backscattering coefficient r (which 
depends on the target atomic number and on the electron primary energy). 

Let us now concentrate our attention on an electron beam impinging on 
a film, i.e. s < R, so that the beam is split into three fractions rjA, Vb 
and rjT representing the fraction of particles absorbed, backscattered and 
transmitted, respectively. Each of them lies in the range [0, 1]. Since the 
total number of particles has to be conserved, then rjA + rjB + t]t = 1 for any 
given thickness. 

Note that if the film is deposited on a substrate, the fractions mentioned 
above are different from those corresponding to unsupported thin films. The 
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difference is due to the backscattering from the substrate. In order to dis- 
tinguish supported and unsupported thin films, we shall use the symbols 
Ca, Cb and Ct for the absorbed, backscattered and transmitted fractions, 
respectively, for supported thin films. Note that each fraction ( lies in the 
range [0, 1]. Conservation of the number of particles for the fractions ( reads 
Ca + Cb + Ct = 1- Obviously, owing to the backscattering from the substrate, 
Cb > Vb and Ct < Vt for any given film thickness. 

The fractions ij and C depend on ^(s), a function related to the scattering 
processes. A rough evaluation of ^(s) is the following: 

C(s) = Natrs , (5.29) 

where N is the number of atoms per unit of volume in the target and atr is 
the elastic transport cross-section. 



Let us observe that 




lim 77^(5) = 0 , 


(5.30) 


Im ?7 b(^) = 0 , 


(5.31) 


lim77T(^) = 1 , 


(5.32) 


lim t]a{0 = 1 - r , 

^ — >-oo 


(5.33) 


lim ysiO = r , 
^—*00 


(5.34) 


lim ?7t(C) = 0 • 


(5.35) 



Let us now consider an electron beam impinging on a film of a given ma- 
terial X of thickness s which has been deposited on a substrate composed of 
a material y (different from material x). Rx indicates the maximum penetra- 
tion range in a bulk target composed of x, and Ry the maximum penetration 
range in a bulk target composed of y. Let us assume that the thickness of 
the substrate is greater than Ry, so that the substrate is a bulk target for 
the primary energy considered. 

The maximum penetration range is a combination of Rx and Ry whose 
value depends on the film thickness. It approaches Ry for s ^ 0 and Rx for 
s > Rx . 

Indicating by the backscattering coefficient of x and by Xy the backscat- 
tering coefficient of y, the following conditions have to be satisfied: 




58 



5 Electrons Impinging on Solid Targets 



lim Ca{0 = 1 - T 

^ — »-oo 

lim Cb{ 0 = rx , 

^ — >-oo 

lim Ct{0 = 0 ■ 

^ — >-oo 



(5.39) 

(5.40) 

(5.41) 



5.2.2 Unsupported thin films 



Let us indicate by an increment of i.e. 

Z\C = ?(s + ^s) - e(s) , (5.42) 

and let us calculate the fraction of particles absorbed by an unsupported 
thin film of thickness s + Z\s. This is given by the fraction absorbed by ^ plus 
the fraction transmitted through ^ and absorbed by A^, plus the fraction 
transmitted through backscattered by A^ and absorbed by and so on, 
in a sequence of infinite reflections: 

7?aK(s + ^s)] =Va{^ + A^) 



= i1a{0 + i1t{0vb{A^)va{0 '^['>1b{0vb{A^)Y 



n—0 



+Vt{0va{A^) ^[t]b{0vb{A0Y 



n—0 



= va(0 



Va{Ovt{Qvb{A^) + rjTiOvAiA^ 
1 - VBiOVBiA^ 



Proceeding in a similar way, we obtain the following: 

Vt(0vb{A^) 



t]b{^ + A^) =r?s(C) 
7 ?t(C + ^0 = 



^-Vb{0vb{A0 ’ 
i1t{0vt{A^) 



^-Vb{0vb{A^) ' 



(5.43) 



(5.44) 



(5.45) 



For 5 ^ oo, riAiO ^ 1 — r, t7b(C) ^ f and ??t(C) 0- Then, as a conse- 

quence, for A^ oo, riA{^ + A^) ^ 1 — r, t]a{A^) ^ 1 — r and t]b{A^) r. 
Let us define 

1 -I- 

fJ- = 



and 



2r 

1 — 
2r 



(5.46) 

(5.47) 
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Using (5.43) for — > oo, it is possible to see that 

riT = ilB- + 1 • 

Let us introduce j3, the derivative of ? 7 s calculated for ^ = 0: 



3 = lim 



Vb{0 



dVBiO 



l«=o ■ 



From (5.44), we obtain the following: 
rjsi^ + A^) - 7 Jb{0 Vb{A0 



vUO 









Therefore 

dVBiO 

d^ 



= (dririO = /3hl(C) - ‘^dVBiO + 1] • 



This equation is equivalent to 

dr]B dr]B ^ n 

TTT\ = 2v(3d^ . 

riB - (1/r) r]B -r 

Remembering that 7ys(0) = 0, the integration of this equation gives 

^ l-exp(-2^/36 

1 - r2 exp(-2z//30 ■ 

As Vt ~ Vb ~ ‘^dVB + 1 and rjA = ^ — riB — ijt-, we can conclude that 



^t(C) = 



1 — 



■exp(-r//3C) 



and 



Va{0 = (!-»') 



1 — r2 exp(— 2i//3^) 

r exp(— 2z//3^) — (1 + r) exp(— z//3^) + 1 



(5.48) 

(5.49) 

(5.50) 

(5.51) 

(5.52) 

(5.53) 

(5.54) 

(5.55) 



1 — r2 exp{—2vP^) 

Equations (5.53), (5.54) and (5.55) have been taken from H.W. Schmidt [4]. 



5.2.3 Supported thin films 

Let us now consider a supported thin film x deposited on a substrate y and 
remember that (^a is the fraction of electrons absorbed in the supported 
thin film, (b the fraction of electrons backscattered and Ct the fraction of 
electrons transmitted across the interface between the film x and substrate 
y. To calculate Ca, let us observe that the fraction rjA is absorbed by the 
film, and then the fraction rjTry comes back across the interface so that the 
fraction rjT'ryrjA is absorbed. From the infinite sum of contributions to the 
absorbed fraction Ca, we obtain 
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Ca =iia + VTfyilA + VrryrjBTyrjA + 

OO 

= r]A l + 

n— 0 
VTTy 



= VA [I + 

V ^-VBTy^ 

In a similar way, we obtain the following, for (^b'- 



CB = r]B + ilT^y '^(ilBry)'^ = — 



Vt^v 



n—0 



1 - VBTy 



We obtain the following for (^t- 



OO OO / -| \ 

^ xn Y^/ '.n — ry) 

Ct = r]T 2 ^[VBry) - r]TTy 2 _^[VBry) = 



(5.56) 



(5.57) 



(5.58) 



n n - VBTy 

n=0 n=0 ^ 

If the material of the film is the same as that of the substrate, we can avoid 
the subscripts x and y in the equations above, so we can write 



Ca = ?7a [^1 + 

C,B = riB + 

Ct = Vt 



VTr 

1 - r]Br 



ri^r 



1 - riBr 
1 — r 



1 - r]Br 



(5.59) 

(5.60) 

(5.61) 



Upon substituting the Schmidt equations (5.53), (5.54) and (5.55) into the 
previous equations, we obtain the following [5]: 



Ca = (1 - r)[l - exp(-i//3C)] , 

Cb = t , 

Ct = (1 - r) exp(-r//3C) . 



(5.62) 

(5.63) 

(5.64) 



Let us observe that the derivative of Ca gives the implantation profile, or 
the depth distribution, of the trapped electrons [5]: 



^ = i^/3(l - r) exp(-z//?C) ^ 

as as 



(5.65) 



Note that the relationship (5.48) between rjT and rjB for unsupported thin 
films of the material x, taking into account the definition of /r (5.46), can also 
be rewritten as 



rjT = 



'(cx - ?7 b)(1 - r^VB) 



(5.66) 
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Substituting (5.66) in (5.57), we can conclude that the fraction of backscat- 
tered particles is given by the following, for a thin film of the material x 
supported on a substrate of the material y [6]: 

CB = VB + {rx-r]B) — -. • (5.67) 

r-x 1 - ryT]B 

All the equations deduced in this section are based on a multiple reflection 
method and were derived neglecting the fact that the absorption, backscatter- 
ing and transmission coefficients change for a given layer owing to the change 
in the energy and angular distribution of the particles that repeatedly are 
incident on the layer. Hence these equations represent an approximate ap- 
proach. However, the equations are certainly useful for gaining insight into 
the subject of the influence of backscattering from the substrate and also to 
understand the theoretical framework. 

As a consequence, the present theory may be informative, but one must 
keep in mind the limits of the analytical procedure. The fact that (5.67) is a 
simple, closed formula certainly renders it attractive. Once the backscattering 
coefficients of the two materials and the fraction of particles backscattered 
by an unsupported thin film are given, the theory allows one to calculate 
the fraction of backscattered particles when the same film is deposited on a 
substrate. 



5.3 Secondary Electrons 

Electron beams impinging on solid targets stimulate the emission of secondary 
electrons, i.e. electrons extracted from the atoms bound in the solid, as a 
consequence of interaction with electrons of the primary beam or with the 
other energetic secondary electrons that travel in the solid. 

Some secondary electrons, following a number of elastic and inelastic col- 
lisions with the atoms of the solid, reach the surface of the solid with enough 
energy to emerge from it. Even some of the primary electrons, after sev- 
eral elastic and inelastic collisions with the target atoms, emerge from the 
surface. So the spectrum of the secondary electrons is contaminated by the 
contribution of the backscattered primary electrons. In the following, we shall 
concentrate our attention on the secondary electrons, choosing to neglect the 
backscattered electrons already treated in the previous sections. In particu- 
lar, we shall describe the theory of the secondary-electron cascade in metals 
as proposed by Wolff in 1954 [7]. 

The process of secondary-electron emission can be conceptually divided 
into two phenomena: the first concerns the production of secondary electrons 
due to interaction between primary electrons and electrons bound in the solid; 
the second is represented by the cascade, whereby the secondary electrons, 
diffusing in the solid, extract new secondary electrons in a cascade process. 
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Since every secondary electron loses energy during its travel within the 
solid, the process proceeds until the secondary-electron energy has decreased 
so that it is not sufficient to extract new secondary electrons or until the 
electron reaches the surface of the solid with sufficient energy to emerge from 
it. 



Let us introduce some relevant physical quantities. We shall indicate by 
Af{r,f2,E,t) the number of electrons at time t between r and r -|- dr, 
between 17 and ^2 + dfl, and between E and E + dE, where 17 represents 
a unit vector in the direction of the velocity of the electron v. X{E) is the 
electron mean free path and E{fl, E; 17', E') the probability that an electron 
at 17', E' is found, after scattering, at 17, E. Let us indicate by S{r, fl, E, t) 
the source, i.e. the density of the secondary electrons due to the bombardment 
of primary particles. The electron cascade is governed by the Boltzmann 
diffusion equation. 



dAf 

dt 



-I- V ■ VAf = 



vAf 



■S 



dE'df2‘ 



\{E') 



E{f2,E; fi',E') (5.68) 



Considering a geometry in which the primary particles are incident nor- 
mally on the target surface, the problem has azimuthal symmetry and in- 
volves the distance to the surface z, the energy E, and the angle between the 
velocity of the secondary electron and the direction normal to the surface, 9. 
If we concentrate our attention on the steady-state condition 



^ = 0 
dt 




(5.69) 


v-VAf = 






+ 


j dE'dn''"'^^''^^;^^'^^'^P{f2,E; f2',E') . 


(5.70) 



Let us indicate by O the angle between 17 and fi' and expand in spherical 
harmonics the three functions Af, E and S\ 



Af{z,cos9, E) = — {21 + l)Af[{z, E)Pi{cos9) , 

47T 

1^0 



(5.71) 



S{z,cos9,E) = —'^{2l + l)Si{z,E)Pi{cos9) , 
^ 1=0 



(5.72) 



E{n,E-f2',E') = P{cos0\E,E') 

- oo 

= —Y,{2l + l)Ei{E,E')Pi{cose) . 
^ 1=0 



(5.73) 
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Let us introduce the function 

= vNi/X{E) , (5.74) 

and note that 

u • VA7 = — y^(2/ + cos0P;(cos0) 

47t ^ oz 






47T 

47T 

47T 



y~^(2/ + 1)^^ cos 9 Pi {cos 9) 

i 

+ l)Ci+i(cos^) + /n_i(cos^)] 

E + A ^ %1(, + DflIeosO) 



dz 



= sE(2'+iWi5) 

i 



( I dipi-i 
\2/ + l dz 



I + 1 d'tpi+i \ 
21+1 dz ) 



P[{cos9) . 



Then, 

M+E)= Si{z,E) + J dE'M+E')Fi{E,E') 



(5.75) 



+ m 



/ I dijji-i 
\2/ + l dz 



I + 1 d-ipi+i \ 
21 + 1 dz ) 



(5.76) 



If we suppose that the secondary electrons are uniformly produced through- 
out the depth, the functions ipi depend only on E and are independent of z. 
If the energy is sufficiently low, the distribution of the secondary electrons 
is spherically symmetric, so one can neglect all the harmonics higher than 
I = 0. Within this approximation, the integro-differential equation (5.76) 
with 5'o = 0 becomes 

p(X) 

ME)= dE'ME')Eo{E,E') . (5.77) 

J E 

Note that, once "00 is known, the secondary-electron energy distribution can 
be calculated as 



j{E) =Nqv = 0oA(P) . 



(5.78) 
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Concerning Fq{E,E'), let us note that 

2tt j sinO d0 E{cos&; E,E')Pk{cos0) 

nTT I 

= 2tt / sinOdO— {21 ~\-l) Fi{E,E') Pi {cos 0)Pk {cos O) 

In 47r 

1^0 

1 

= -j2i2i+m{E,E') 

^ 1=0 

1 °° s> 

1=0 

= Ek{E,E'), (5.79) 

so that 

Eo{E,E') =27ry sin 6 > J^(cos 0; ii') . (5.80) 

Therefore, within the approximation of spherical symmetry of the secondary- 
electron distribution and ignoring all the harmonics higher than I = 0, (5.77) 
can also be written as the following: 

poo pTT 

V'o(-E) = 27 t / dE''iljo{E') sin0d0F{cos0;E,E') . (5.81) 

J E Jo 

Here Fq{E,E') represents the total probability of scattering between the 
energies E and E', independently of the angle. 

Let us assume that, on average, electrons of energy E' < 100 eV lose 
half of their energy at each collision. Let us introduce the average energy 
after scattering, E, defined by 

E = -fE' . (5.82) 

Here E' is the energy before the scattering and, for energies four times greater 
than the Fermi energy, 7 is ~ 1/2. Within this approximation, 

Fo(F;, E') =2S{E - E) =2 6{E - -fE') . (5.83) 



Pi{u)Pk{u) du 



Note the factor of 2, which takes into account the fact that there are two 
electrons after each electron collides. 

Equation (5.77) therefore becomes 



il^oiE) = 2js{E- jE')ME') dE' , 



(5.84) 



which, as the reader can easily verify, produces the following solution: 



'il)o{E) 




E 

pE). 



(5.85) 
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Let us now define the function x{E), so that 

= 1 . (5.86) 

Our solution is proportional to E~^: 

ME) oc E~^^^'> . (5.87) 



Assuming that the validity of the Wolff theory can be extended up to 
~1000 eV, the semi-empirical law proposed in 1977 by Sickafus [8, 9, 10] 
immediately follows. As the inelastic mean free path X{E) is approximately 
proportional to E/logE, i.e. 



A oc 



E 

log if 



(5.88) 



we can see that the secondary-electron energy distribution is approximately 
given by 



C C E~'^ 

j{E) = il}oKE) = ^ logE ’ 

where C is a constant and m = x — 1. Then, 
log j{E) = log C — log(log if) — m log if . 



(5.89) 

(5.90) 



Note that the function log(logif) is almost constant. As a result, in the range 
of energies 10-1000 eV above the low-energy cascade peak. 



j-(if) = Aif— , 



(5.91) 



where A is almost constant and m ~ 1. 



5.4 Comparing Theory and Experimental Data 

We have introduced many physical quantities and described them through 
semi-empirical approaches in this chapter. Today, the most accurate theo- 
retical calculation of these physical quantities can be performed by Monte 
Carlo simulations, which will be faced in the next chapter. Consequently, the 
present chapter shows only a comparison of the results obtained from the 
multiple reflection method with the experimental data available concerning 
backscattering from supported thin films. This has been done because, despite 
the simplicity of the approach, there is substantial agreement between theory 
and experiment. The other quantities introduced in this chapter (backscat- 
tering coefficient, depth distribution of the trapped electrons and positrons, 
secondary-electron emission, and so on) will be calculated in the next chap- 
ter by introducing the elastic and inelastic cross-sections (described in the 
previous chapters) into the Monte Carlo simulations. 

In order to investigate the applicability of the multiple reflection method 
and to compare its results with the experimental data reported by Niedrig in 
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Fig. 5.1. Electron backscattering ratio ("s of Ag surface films deposited on A1 
bulk substrates versus surface film thickness z. The data presented concern 20 keV 
electron beams irradiating the targets in the +z direction. +: theoretical calculation 
based on the multiple reflection method (5.67) [6]. •: Niedrig’s experimental data 
[ 11 ] 




.(A) 



Fig. 5.2. Electron backscattering ratio Cs of An surface films deposited on A1 
bulk substrates versus surface film thickness The data presented concern 20 keV 
electron beams irradiating the targets in the +z direction. +: theoretical calculation 
based on the multiple reflection method (5.67) [6]. •: Niedrig’s experimental data 
[ 11 ] 
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Fig. 5.3. Electron backscattering ratio of An surface films deposited on A1 
bulk substrates versus surface film thickness z. The data presented concern 30 keV 
electron beams irradiating the targets in the -\-z direction. +: theoretical calculation 
based on the multiple reflection method (5.67) [6]. •: Niedrig’s experimental data 
[ 11 ] 




[11] we used Niedrig’s experimental data for the values of fsiz), D and rg in 
(5.67). The agreement between the theoretical and the experimental values 
of (b{z) is excellent (see Figs. 5. 1-5. 3). 

Note that (5.67) was derived by the multiple reflection method, assum- 
ing that fixed values for the probabilities of absorption, backscattering and 
transmission could be attributed to both the surface layer and the substrate. 
However, this is not strictly valid, as these probabilities depend on the an- 
gular distribution of the particles and on their energies. Since the particles 
dissipate energy in ionisation, electron excitation, plasmon emission and so 
on, the assumption of fixed probabilities represents an approximation [12]: 
our comparisons were motivated by the necessity to evaluate the accuracy of 
such an approximation. 

Despite the approximations involved, the agreement, both in the trend 
and in the absolute values, with Niedrig’s experimental data concerning the 
values of Cb{z) shows that the multiple reflection method is a reliable ap- 
proach for calculations of backscattering from supported thin films and that 
it is useful for fast and quite accurate evaluations when data concerning un- 
supported thin films and backscattering coefficients are available. 

In conclusion, the accord that we found indicates that (5.67) can be safely 
used for calculating backscattering from supported films and that the multiple 
reflection method is rather a good approximation. 
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6 Monte Carlo Simulations 



6.1 The Monte Carlo Method 

The study of the physical quantities that are involved in electron-solid 
interaction (backscattering coefficient, energy and angular distribution of 
backscattered electrons, electron penetration in thin solid films, electron im- 
plantation profiles in bulk targets, secondary-electron energy distribution, 
etc.) requires approximations and approaches which are valid only for lim- 
ited ranges of energy and angle and for selected materials. All the analytical 
approaches we have proposed involve approximations and have been deduced 
using semi-empirical models. 

A completely different approach to the determination of the quantities 
mentioned above is represented by the Monte Carlo method. Within statis- 
tical uncertainty, which can be reduced arbitrarily, this method of simulat- 
ing electron-solid interactions is very accurate and gives results in excellent 
agreement with the experimental data available. 

The Monte Carlo method is a mathematical tool utilising random num- 
bers. This method can obtain accurate results about the absorption, backscat- 
tering and transmission of electrons penetrating both supported and unsup- 
ported thin films, and can also provide implantation profiles of the absorbed 
electrons, the energy and angular distributions of backscattered and trans- 
mitted electrons, and spectra of secondary electrons. The cross-sections and 
probabilistic laws for interaction of electrons with the atoms constituting the 
target are known, so that they can be used in the Monte Carlo code. 

When a stream of electrons impinges on a solid target, the incident par- 
ticles (during their travel within the solid) lose energy and change direction 
at each collision with the atoms bound in the solid. Nuclear collisions deflect 
the electrons without relevant kinetic-energy transfer, owing to the large dif- 
ference between the masses of the electron and the nucleus: the differential 
elastic scattering cross-section has to be used in order to describe these pro- 
cesses (see Chap. 3). 

On the other hand, excitations and ejections of atomic electrons, and 
plasmon excitations affect the energy dissipation but only slightly affect the 
direction of the electron in the solid. These can be described as inelastic 
processes that are essentially governed by the equations for the stopping 
power and/or the inelastic cross-section (see Chap. 4). 
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All the cross-sections, mean free paths and stopping powers can thus be 
accurately calculated by the Monte Carlo method in order to obtain the 
macroscopic characteristics of the interaction processes (see Chap. 5), by 
simulating a large number of single particle trajectories and then averaging 
them. 

However, as the Monte Carlo method is statistical, the accuracy of its 
results depends on the number of simulated trajectories. Recent evolution in 
computer calculation capability means we are now able to obtain statistically 
significant results in very short times of calculation. 



6.2 Random Variables 



We need to introduce the concept of a random variable in order to understand 
the Monte Carlo method. Let us begin with the discrete random variable. In 
order to specify a discrete random variable one has to provide the n values 
that ^ can assume, with their probabilities of occurrence. 

In other words, we have to specify and their probabilities 

Pi = P{i = Ci))P 2 = p{^ = 6)) ■■■,Pn = p{^ = Cn)- Here, we have indicated 
by the probability pi that a value of ^ chosen at random is equal to 

^i, with i = 1, 2, . . . , n. Of course, X)r=i Pi ~ for every z = 1, 2, . . . , n. 

Pi > 0. Once a discrete random variable has been specified, its expected value 
can be calculated from 

= = ( 6 . 1 ) 

E^p^ V 



In order to define a continuous random variable one has to give the 
range of variability and the distribution of the probabilities. If a continuous 
random variable ^ is defined, this means that we have specified a real interval 
(a, b) and a probability density p{x) of the distribution of the variable. Let 
us now consider the real interval (c, d) in (a, b). If P{c < ^ < d) indicates the 
probability that a value of ^ taken at random falls in the interval (c, d ) , then 



P(c < ^ < d) = f p{x) dx . (6.2) 

J c 

As fl^pix) dx = 1, the expected value of ^ can be calculated from 

(^) = f xp{x) dx . (6.3) 

J a 

If we need to calculate the expected value of a function /(^), this is given by 
ifiO) = [ f{x)p{x) dx . 

J a 



(6.4) 
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6.2.1 Random Variable Uniformly Distributed 
in the Interval (0, 1) 



A random variable fj, is said to be uniformly distributed in the interval (0, 1) 
if it is defined in the range (0, 1) and if its probability density is simply given 
by Pf,{x)=l. 

The probability that a random variable uniformly distributed in the in- 
terval (0, 1) assumes a value between c and d (where c > 0 and d < 1) is the 
length of the interval. In fact, 

pd nd 

P{c < p < d) = / Pfi{x)dx= / dx = d — c . (6.5) 

J c j c 

The expected value of a random variable uniformly distributed in the interval 
(0, 1) equals 1/2, as follows: 



(/r) = / xpfj,{x) dx = / xdx 
Jo Jo 



( 6 . 6 ) 



6.2.2 Random Variable Distributed in a Given Interval 
with a Given Probability 

If f indicates a random variable defined in the interval (a, b) and distributed 
with a given probability density p(x), and p a random variable uniformly 
distributed in the range (0, 1), the values of ^ are related to those of p and 
can be obtained from the equation 

f p{x) dx = p . (6.7) 

J a 

Let us consider the function /(x) defined by 



/(x) = f p{z) dz . 
J a 



( 6 . 8 ) 



The function /(x) has the following properties: it increases when x in- 
creases because /(a) = 0 and f{h) = 1 and, by the definition of probability, 
df{x)/dx = p{x) > 0. 

To proceed, let us consider the straight line 

y = p, (6.9) 

where p is & real number in the interval (0, 1). The line is parallel to the 
X axis and must intersect the increasing function /(x) in only one point. 
Therefore, a intersection between the two curves exists and occurs only once, 
i.e. for each value of p, (6.7) has a solution, and this solution is unique. 

Let us now show that the distribution f defined as described here has a 
probability distribution p{x). If we consider an interval (c, d) contained in 
the interval (a, b), then 

P(c < 5 < d) = P[/(c) <p< f{d)] . 



( 6 . 10 ) 
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As the probability density of ^ equals 1, then 

p[fic)<^i<f{d)] = m-fic) ( 6 . 11 ) 

and, on the other hand, 

f{d)-f{c)= ( p{x)dx. (6.12) 

J C 

From the last three equations, we obtain the result that 

P{c < ^ < d) = j p{x) dx , (6.13) 

and hence the random variable ^ is distributed with a probability density 
p{x) in the interval (a, 6). 



6.2.3 Random Variable Uniformly Distributed 
in the Interval (a, b) 



We are now interested in calculating a random variable rj uniformly dis- 
tributed in the real interval (a, b), which is useful in Monte Carlo simula- 
tions of electron-solid interactions. The requirement that the distribution is 
uniform implies that the corresponding probability density is given by 

Pvi^) = ) (6-14) 



so that if we indicate by /i a random variable uniformly distributed in the 
range (0, 1), then rj must satisfy the following equation: 




As a consequence, rj = a + p{b — a) and its expected value is given by (77) = 
(a -|- &)/2. 



6.2.4 Random Variable with Poisson Distribution 

The Poisson distribution is another very important distribution for Monte 
Carlo simulations of electron-beam interactions with solid targets. This dis- 
tribution is defined by the following probability density: 

= ^exp , (6.16) 

where A equals a constant (whose meaning will be explained later). A random 
variable x defined in the interval (0,oo), whose probability density is the 
Poisson distribution, is obtained by the solution of the equation 




(6.17) 
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where /i indicates a random variable uniformly distributed in the range (0, 1). 
As a consequence, we can conclude that 

X = -Aln(l - Ai) ■ (6.18) 

As the distribution of 1 — /i equals that of /i, we can write 



X = -Aln(/x). (6.19) 

It is therefore possible to see that the constant A is the expected value of y: 

(x) = J dx = j j exp dx 

= [-xexp(-x/A)]|)f’ + / exp(-cc/A)i 
Jo 



I dx 



= — Aexp(— x/A)|jf’ = A . 



( 6 . 20 ) 



6.2.5 Pseudo- Random-Number Generators 

If we are able to generate random numbers uniformly distributed in the range 
(0, 1), it is possible, as we have previously shown, to generate random num- 
bers distributed with any given probability density in any given interval. 
In other words, we can say that the uniformly distributed random variable 
underlies the generation of any other random number. As people generally 
use computer procedures that generate sequences of pseudo-random num- 
bers uniformly distributed in the interval (0, 1), the problem arises of check- 
ing the uniformity of a given pseudo-random-number generator. A simple 
way to check the uniformity of a given distribution is to simulate the value 
of 7T = 3.14 .. . . This can be performed by generating numerous pairs of 
pseudo-random numbers distributed in the range (—1, 1). If /i is a sequence 
of pseudo-random numbers distributed in the interval (0, 1) then 2/r — 1 is 
a sequence of pseudo-random numbers distributed in the interval (—1, 1). 
If the number of pairs is statistically significant and the distribution is uni- 
form, then the ratio between the number of pairs that lie in the unit circle 
and the total number of generated pairs should approach the value 7t/4. In 
other words, a comparison between the Monte Carlo-calculated area of the 
unit circle and its exact value can be used as a method to check the level of 
uniformity of a pseudo-random- number generator. 



6.3 A Simple Monte Carlo Scheme 

We can introduce a simple, reliable procedure, whereby a flux of electrons 
with energy Eq impinges on a homogeneous target in the +z direction. The 
electrons of the beam are elastically scattered by the atoms of the target, and 
lose their energy in inelastic collisions with the atomic electrons and through 
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plasmon excitations. For the simple procedure we are now describing, we 
assume that all the energy lost by the electrons of the incident beam can 
be described by utilising a continuous-energy-loss approximation. Electrons 
are thus assumed to lose energy continuously inside the solid, and the energy 
loss processes are all incorporated into the stopping power. The statistical 
fluctuations of the energy loss are therefore disregarded in this case. 

The elastic scattering cross-sections are calculated assuming that the first 
Born approximation is valid throughout the electron travel in the solid target. 
The systematic errors introduced by this approximation are expected to be 
more important for the lowest primary energies. 

We adopt spherical coordinates (r, 0, (j)), and the path-length distribution 
is assumed to follow Poisson statistics, with a mean value equal to the electron 
elastic mean free path. The step length As is then given by 

Z\s = -Aei In(^i) , (6.21) 

where fXi is a random number uniformly distributed in the range (0, 1) while 
Xei represents the elastic mean free path. This is given by 

where N is the number of atoms per unit of volume in the target and aei is 
the total elastic scattering cross-section 

CTei = 27 t / -—sinOdO. (6.23) 

Jo ds2 

For a Wentzel-like potential V{r) = — (Ze^/r) exp(— r/a), in the first Born 
approximation, the total elastic scattering cross-section is given by 



TrZ'^e'^ 1 


(6.24) 


“ E2 a(2 + a) ’ 


so that the elastic mean free path is 




a(2 -I- a)E^ 

~ Nne^Z^ ’ 


(6.25) 


where 






(6.26) 





The polar scattering angle 0 after an elastic collision is selected assuming 
that the probability P{6) of elastic scattering into an angular range from 0 to 
0 is a random number /i 2 , uniformly distributed in the range (0, 1). In other 
words, the integrated probability for scattering in the angular range from 0 
to 9 equals fj, 2 - 

2tt J^(da/dn) sinidd'd 
2tt J^{da/df2) siniddiJ 



d2 = P{9) 



(6.27) 
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From the equation 



P{9) = 


(1 -1- a/2)(l — COS0) 


(6.28) 


1 — cos 9 + a 


it follows that 




cos 9 = 


2a/X2 


(6.29) 




2 + 0; — 2/7,2 



The azimuthal angle </> can take on any value in the range (0, 27 t) and is 
selected by a random number rj uniformly distributed in the interval (0, 27 t). 

Note that both the polar scattering angle 6 and the azimuthal angle (j) 
are calculated relative to the last direction in which the electron was moving 
before the collision. The angle 6'^ at which the electron is moving after the 
last deflection, relative to the z direction, is given by 

cos 9'^ = cos 9z cos 9 + sin 9^ sin 9 cos (j) , (6.30) 



where 9z is the angle relative to the 2 direction before impact. The motion 
Az along the z direction can then be obtained from 

Az= As cos 9'^. (6.31) 



The new angle 9'^ is the incident angle 9^ for the next path length. Electrons 
are followed into the solid target until their energy becomes lower than the 
mean ionisation energy J (in eV, J = 13. 6Z for Z < 10; otherwise J = 
(9.76 + 58.8Z”^-^®)Z) or until they emerge from the surface (backscattered 
particles) or from the back of the target (transmitted particles). In such a way, 
we avoid taking into account very low energies, for which the elastic-scattering 
calculation requires the partial-wave expansion method. This procedure is 
quite accurate if the primary energy of the impinging electron beam is higher 
than ~10 keV. 

For these primary energies, the energy loss can be safely calculated using 
the Bethe formula [1], 



dE 

ds 



2ne^NZ , /1.166E\ 



E 



■ In 



J 



(6.32) 



or with semi-empirical expressions, such as the following, 

dE K^NZ^/^ 

~ ~dl~ E2/3 ’ 



(6.33) 



proposed in 1972 by Kanaya and Okayama (with K^. = 360 eV®/^ A^) [2]. 
This last formula allows one to analytically evaluate the maximum range of 
penetration, whereby 



R = 




dE 

dE/ds 



3Ep/^ 

5KeNZ^/^ ■ 



(6.34) 
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Once the segment of the trajectory As has been calculated, the energy 
loss AE can be approximated by 



.L as 



dE 

ds 



As . 



(6.35) 



Note that, with this approach, we neglect the energy loss straggling, i.e. 
the statistical fluctuations of the energy losses. 

So, the simple Monte Carlo simulation described here can be safely utilised 
for primary energies higher than 10 keV in the calculation of the backscatter- 
ing coefficient, for example [3, 4]: the backscattering coefficient is calculated 
by counting the number of electrons that emerge from the surface of a bulk 
target and dividing it by the total number of simulated trajectories (which 
should be more than 10 000 in order to obtain statistically significant values 
of the backscattering coefficient). 



6.4 A More Sophisticated Simulation 



The Monte Carlo scheme previously described can also be utilised for energies 
lower than 10 keV. To do this, it is necessary that the cross-sections utilised to 
describe the elastic scattering are calculated with the partial-wave expansion 
method, and the stopping power and the inelastic mean free path are com- 
puted utilising optical data and the dielectric function. With the previously 
described scheme, it is possible to simulate the interaction of electrons of low 
to medium primary energy (1-10 keV) and to extend each trajectory until 
the electron energy becomes lower than ~50 eV, with these cross-sections 
and stopping powers. 

The following analytical expression can be utilised, concerning the use of 
the relativistic partial-wave expansion method, to expedite the calculations 
and to approximate the differential elastic scattering cross-section [5, 6, 7]: 



da 

dfi 



<P{Z, E) 



1 

(1 -cos6»-HT)2 • 



(6.36) 



Note that, with <!>{Z,E) = Z^e'*/4if2 and T{Z,E) = 0, this equation be- 
comes identical to the classical Rutherford formula (3.38), while with <P{Z, E) 
= Z^e'*/4i?2 and T(2’, £1) = a{Z,E) = (me^7r^//i^)(Z^/^/if), it becomes 
the screened Rutherford formula (corresponding to a Wentzel-like potential) 
(3.36). So, as (6.36) is identical to the equation we used in the previously 
described Monte Carlo scheme, we can use the same formal structure as 
in that procedure. However, it is necessary that <P{Z, E) and T{Z, E) are 
computed in order to obtain, by using (6.36), the same values of the total 
cross-section aei and of the transport cross-section atr as calculated by util- 
ising the quantum-relativistic partial-wave expansion method. [8, 9, 10, 11]. 
Note that these quantities are defined, respectively, by the following: 
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/ d(j 
(1 — cos 6)— —dfl . 
all 



(6.37) 

(6.38) 



Therefore, (6.36) allows one to sample the polar scattering angle with a 
closed formula similar to that used in the previous Monte Carlo scheme: 

2TP{9) 



cos 0=1 — 



(6.39) 



2 + r-2P(0) ■ 

Let us now show how to calculate d>{Z,E) and T{Z,E) by use of our 
knowledge of the total elastic scattering cross-section agi and of the transport 
elastic scattering cross-section atr mentioned above. Note that from (6.37) it 
follows that 

, (6.40) 

so that the differential elastic scattering cross-section can be rewritten as 

^_r(r + 2) 

dfi 



(6.41) 



47t (1 — COS0 -I- TY 
Using (6.38) and (6.41), we obtain, concerning the transport cross-section 



(j^r — (7 el 



T{T + 2) fr + 2 



In 



r 



- r 



(6.42) 



Let us calculate now the ratio n between the transport elastic scattering 
cross-section and the total elastic scattering cross-section: 



^ tr 
^ el 



T + 2 



• In 



r 4-2 
T 




(6.43) 



The values of atr and of aei have been numerically calculated, so that it is 
possible to find S’ as a function of Z and E. The values of T as a function 
of Z and of the electron energy E can be subsequently computed, by use of 
our knowledge of S, by utilising a bisection algorithm. Once the values of 
T{Z,E) have been stored in a file, it is possible to calculate the value of T 
corresponding to the particle energy at every step of the electron trajectory 
during the simulation. The sampling of the scattering angle 9 can be easily 
performed by inserting the value of T (found using an interpolation algorithm 
such as a cubic spline) into (6.39) and selecting a random number uniformly 
distributed in the range (0, 1) for P{9). 



6.4.1 Surface Films 

The interface between the film and the substrate must be taken into account 
for surface films. After crossing the interface, the change in the scattering 
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probabilities per unit length in passing from the film to the substrate and 
vice versa has to be considered. Let us denote by pi and p 2 the scattering 
probabilities per unit length for the two materials, where pi refers to the 
material in which the last elastic collision occurred and p 2 to the other mate- 
rial. d indicates the distance along the scattering direction between the initial 
scattering event and the interface. If rnd is a random number uniformly dis- 
tributed in the range (0, 1), the step-length Al is given by 



(l/pi)[-ln(l - rnd)] 
for 0 < rnd < 1 — exp(— pid) 
d + {\/p 2 )[- ln(l - rnd) - pid] 
for 1 — exp(— pid) < rnd < 1 



(6.44) 



If Pi = p 2 (i.e. if there is no interface), the interface is only an imaginary 
boundary at a depth z below the surface and, for each value of rnd, Al = 
Ai[— ln(l — rnd)]. If d is very large with respect to 1/pi, then 1 — exp(— pid) « 
1 and, as a consequence, Z\/ « (l/pi)[— ln(l — rnd)] = Ai[—ln(l — rnd)]. In 
other words, when the particle is far away from the interface, its behaviour 
is practically the same as that described by (6.21). Note that, in general, if 
rnd < 1 — exp(— pid) then the next elastic collision will occur in the same 
material and the interface will not be crossed. On the other hand, when 
rnd > 1 — exp(— pid), then the next elastic collision will be in the other 
material because the interface will be crossed. In this case, the particle will 
travel a distance d along the scattering direction in the material where the 
last collision has occurred according to the scattering probability for that 
material and, after crossing the interface, it will travel a distance Al — d 
along the scattering direction according to the scattering probability of the 
other material. The energy loss is approximated by utilising the same random 
number rnd used to calculate Al: 



{dE/dl)i{l/pi)[— ln{l — rnd)] 
for 0 < rnd < 1 — exp(— pid) 

{dE/dl)id + {dE / dl) 2 (l/p 2 )[— ln(l — rnd) — pid] 
for 1 — exp(— pid) < rnd < 1 



(6.45) 



where {—dE/dl)i is the stopping power in the material where the last elastic 
collision occurred and {—dE/dl )2 the stopping power in the other material. 
Note that this approach is correct only if the film thickness is greater than 
the electron elastic mean free path, which means that the (6.44) and (6.45) 
cannot be used when the film thickness is smaller than the electron elastic 
mean free path. 



6.5 Another Monte Carlo Scheme 

We shall now briefly describe another Monte Carlo scheme based on the 
classical-trajectory approximation. This scheme distinguishes the elastic and 
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the inelastic collisions during electron travel inside the solid. The step length 
As is given by 

As = -At ln(/xi) , (6.46) 

where fXi is a random number uniformly distributed in the range (0, 1) and 
At is the total mean free path: 

A. = ■ _r ■ _r . (6.47) 

\nel '^el 

Before each collision, a random number uniformly distributed in the range 
(0, 1) is generated and compared with the probability of inelastic scattering 
Pinei- The probability of inelastic scattering is given by 

Pinel = \nel! ^\nel ^el ) ’ (6.48) 

while, of course, that of elastic scattering is Pei = 1 — Pinei- If the random 
number generated is less than or equal to the probability of inelastic scatter- 
ing, then the collision will be inelastic; otherwise, it will be elastic. 

If the collision is elastic, then the calculation of the polar scattering angle 
proceeds as in the previously described Monte Carlo scheme. The azimuthal 
angle is calculated as a random number uniformly distributed in the range 
(0, 2tt). 

If, on the other hand, the collision is inelastic, then the energy loss is 
computed by utilising a random number p 2 , uniformly distributed in the 
range (0, 1), so that 

Jo[dX-!jd{AE')]d{AE') ’ 

where AE is the energy transferred to a secondary electron. The secondary 
cascade may then be followed. 

6.5.1 Angular Deflection in Electron Electron Collisions 

Note that, in the elastic electron-electron collisions corresponding to inelastic 
electron-atom scattering, the electrons not only lose energy but also suffer 
angular deflection. 

To evaluate the angular deflections suffered by the electrons when the 
electron-atom collision is inelastic, we shall introduce the classical binary- 
collision model, which is sufficiently accurate for many practical purposes. In 
other words, we shall treat the scattering problem of the angular deflection in 
an inelastic electron-atom collision by studying an elastic electron-electron 
classical collision. 

Let us consider an electron travelling along the ^ direction, and suppose 
that a free electron is at rest at a distance r. Let us indicate by vq the velocity 
of the incident electron (the primary electron) before impact, by v its velocity 
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after impact, by u the velocity after impact of the electron initially at rest 
(the secondary electron), and by 9 and "9 the scattering angles of the primary 
and of the secondary electron, respectively, after impact. 

Energy and momentum conservation imply that 



v1 = , (6.50) 

?;o = "c cos 0 + u cos d , (6.51) 

rising — t6 sin -d = 0 . (6.52) 

Therefore, 

cos^ 9 = {vq — ucos'dY , (6.53) 

sin^ d = sin^ -d (6.54) 

and, as a consequence, 

= 2uwo cos d . (6.55) 

Then, 

u = ?;ocosd , (6.56) 

= t,2(i_cos^d) = vgsin^d . (6.57) 



Moreover, 



sin^ 9 = — sin^ d = 



Vq cos^ d 
Wg sin^ d 



sin^ d 



so that 



(6.58) 



sin^ 9 = cos^ d . 



(6.59) 



The energy of the primary electron is E = mv^j2. Let us indicate by AE the 
energy lost by the primary electron and transferred to the secondary electron: 



AE = 



1 2 1 2 
- 2 ™^ 



2 

-mw . 



(6.60) 



Straightforward algebraic manipulations allow us to conclude that 

AE 2 a •2/1 
— — = COS d = sin 9 . 

hj 



(6.61) 



6.5.2 Secondary Electrons 

The incorporation of secondary electrons requires a knowledge of their ini- 
tial angles and energies. If 9 and </> are respectively the scattering and the 
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azimuthal angle of the primary electron, using the classical binary-collision 
theory, the scattering angle t? of the secondary electron should satisfy the 
relationship 

cos-d = sin0, (6.62) 

and the azimuthal angle (p of the secondary electron is given by 

ip = TT + 4> . (6.63) 

The last two equations are used to represent the ionisation processes. 
Since slow secondary electrons are generated with spherical symmetry, in the 
case of the Fermi sea excitation we can assume a random direction of the 
secondary electrons: 

'& = TTP3 , (6.64) 

p = 27T/i4 . (6.65) 



Here ps and p 4 are random numbers uniformly distributed in the range (0, 1). 
The primary-particle energy loss AE is transferred to the secondary electron. 
The electron that has less energy is generally considered as the secondary. 

The surface energy barrier influences the energy distribution of the low- 
energy electrons, so that secondary electrons cannot emerge from the surface 
if they have an angle 0 with respect to the normal higher than 



= cos 



-1 



I Ef + U 
E + Ef + U 



( 6 . 66 ) 



where Ef is the Fermi energy and U the work function. 



6.6 Comparing Theory and Experimental Data 

In the previous chapters, we showed that, for scattering angles higher than 
5°, the accuracy of the calculation of the differential elastic scattering cross- 
sections presented was 1-2%, while that for the total elastic scattering cross- 
section was of the order of 5-6%. On the other hand, the total cross-section 
does not feature strongly in multiple-scattering processes: a systematic study 
of the first transport elastic scattering cross-section has shown that the ac- 
curacy of the present approach is of the order of 2%. Taking into account the 
inaccuracies in (i) the evaluation of the elastic and inelastic cross-sections, 
(ii) the approximations introduced by using cubic spline interpolations of 
tabulated data, (iii) the stopping power used to calculate the energy losses 
and (iv) the statistical uncertainty of the Monte Carlo procedure, we are 
confident in concluding that our simulations give results with an accuracy of 
5-15% for electron and positron energies higher than ~100 eV. In the case 
of lower energies, we may anticipate greater inaccuracies because the evalu- 
ation of the cross-sections is less accurate. In particular, in this low-energy 
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regime (energies lower than 100-200 eV), the accuracy of the simulation 
should be evaluated by directly comparing the results with the experimental 
data available, as it is well known that theoretical evaluation of the relevant 
cross-sections is still an open problem. 

In Fig. 6.1, we compare the Monte Carlo calculation results for the frac- 
tions of electrons absorbed by and transmitted through unsupported thin 
films of Cu of different thicknesses with the experimental data of Cosslett 
and Thomas [12]. In Figs. 6.2, 6.3, the electron backscattering coefficients for 
5 keV and 10 keV electrons, respectively, impinging on various targets, are 
presented and compared with experimental data. The agreement between the 
Monte Carlo results and the experimental data considered here is satisfactory, 
with an accuracy of the order of 1-10%. 

The depth distribution of 3 keV electrons absorbed in a 100 A thin film of 
gold deposited on a silicon dioxide substrate is shown in Fig. 6.4, as calculated 
with the Monte Carlo code. 

The fraction of electrons backscattered by a supported thin film Cb{z) as 
a function of the film thickness z can be directly calculated by a Monte Carlo 
simulation. On the other hand, by using Monte Carlo data concerning the 
fraction for unsupported films of thickness z and the backscattering 

coefficients rt and of the materials constituting the film and the substrate, 
(5.67) can also be used to calculate Cb{z). A comparison between the values 
of (b (z) obtained by these two different procedures allows us to further check 




z(A) 



Fig . 6 . 1 . Fractions t]a and rjT of electrons absorbed and transmitted, respectively, 
for unsupported thin films of Cu as a function of the film thickness z. The primary 
energy is 10 keV. •: rjA, Cosslett and Thomas experimental data [12]. -|-: rjA, Monte 
Carlo data [13]. o: rix, Cosslett and Thomas experimental data [12]. x: t^t, Monte 
Carlo data [13] 
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Fig. 6.2. Backscattering coefficient r as a function of the atomic number Z for 
streams of 5 keV electrons irradiating targets in the +« direction. +: Monte Carlo 
data [13]. •: Bishop experimental data [14]. o; Hunger and Kiichler experimental 
data [15] 
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Fig. 6.3. Backscattering coefficient r as a function of the atomic number Z for 
streams of 10 keV electrons irradiating targets in the +z direction. +: Monte Carlo 
data [13]. •: Bishop experimental data [14]. o; Hunger and Kiichler experimental 
data [15] 




84 



6 Monte Carlo Simulations 




.(A) 



Fig. 6.4. Depth profile P{z) of 3 keV electrons penetrating into a thin film of An 
deposited on a Si02 substrate. Integration of P{z) from z = 0 to z = R gives the 
absorption coefficient 1 — r. The An film thickness is 100 A. o: Monte Carlo data 



the accuracy of the multiple reflection method for studying electron backscat- 
tering from supported thin Aims. Such a comparison is shown in Figs. 6.5 and 
6.6, where we show the result of simulating 5 and 10 keV electron penetration 
in gold thin Aims deposited on carbon substrates. The agreement is satisfac- 
tory: the differences between the direct Monte Carlo computations and the 
results obtained using (5.67) are always smaller than ~10-15% of the Monte 
Carlo data. Also, the general trend as a function of the thickness is what one 
expects on the basis of physical considerations. 

Note that inelastic effects are treated phenomenologically in the Monte 
Carlo approach: they are included in the analytical approach with exactly 
the same level of accuracy as in the Monte Carlo because the values of 
rt and utilised in (5.67) are the same as those used in the code. However, 
even though we have used values of t]b{z), rt and Vg calculated by the Monte 
Carlo method in (5.67), note that the directly calculated Monte Carlo values 
of Cb{z) were obtained through independent computations. 

We present our results concerning the stopping profiles of electrons im- 
planted in silicon dioxide and of positrons implanted in copper, respectively, 
in Figs. 6.7 and 6.8. The primary energies considered were 3.0, 5.0 and 10.0 
keV. The stopping profiles have been given as RP{z,Eo), where R is the 
range and 2 is the depth inside the solid measured from the surface. Note 
that, for any given, fixed primary energy Eg, integration of P{z,Eq) from 
z = 0 to z = R gives the absorption coefficient. Each curve of implanta- 
tion was obtained by simulating 10® electron or positron trajectories. We 
also performed similar calculations for carbon, silver and gold. Note that it 
is possible to calculate the backscattering coefficient r and the range R as 
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Fig. 6.5. Electron backscattering ratio of An surface films deposited on C bnlk 
substrates vs surface film thickiress 2 . The data presented concern 5 keV electron 
beams irradiating targets in the +z direction. +: theoretical calculation based on 
the multiple reflection method (5.67) [13]. •: Monte Carlo data [13] 




2(A) 



Fig. 6.6. Electron backscattering ratio of An surface films deposited on C bulk 
substrates vs surface film thickness 2 . The data presented concern 10 keV electron 
beams irradiating targets in the +z direction. +: theoretical calculation based on 
the multiple reflection method (5.67) [13]. •: Monte Carlo data [13] 
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Fig. 6.7. Stopping profiles RP{z,Eq) vs z (where R is the range, P{z,Eo) the 
depth distribution, « the depth inside the solid measured from the surface and Eq 
the primary energy) for electrons in Si02- For any given primary energy, integration 
of P{z, Eo) from z = 0 to z = R gives the absorption coefficient 1 — r. The data 
presented concern Monte Carlo simulations of electron beams irradiating targets in 
the +z direction [7]. +: 3 keV. •: 5 keV. o: 10 keV 
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Fig. 6.8. Stopping profiles RP{z,Eq) vs z (where R is the range, P{z,Eq) the 
depth distribution, a the depth inside the solid measured from the surface and Eq 
the primary energy) for positrons in Cu. For any given primary energy, integration 
of P{z, Eo) from z = 0 to z = R gives the absorption coefficient 1 — r. The data 
presented concern Monte Carlo simulations of positron beams irradiating targets 
in the +z direction [16]. +: 3 keV. •: 5 keV. o: 10 keV 
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a function of the primary energy and of the atomic number, starting from 
the implantation profiles. The results reported compare excellently with both 
calculated and experimentally obtained data of other investigators. With re- 
gard to positrons, the calculated stopping profiles are an essential first step in 
solving the transport equation that is used in the application of the positron 
annihilation technique. The collisional processes of positron beams impinging 
on solids have recently received great attention, owing to the possibility of 
performing non-destructive investigations of point and extended defects of 
surfaces, interfaces and bulk materials by positron annihilation spectroscopy. 

In conclusion, let us consider the problem of the secondary-electron emis- 
sion from solids irradiated by a particle beam. It is an important problem, 
particularly in connection with analytical techniques which utilise secondary 
electrons to investigate chemical and compositional properties of solids in the 
near-surface layers, namely Auger electron spectroscopy and X-ray photo- 
electron spectroscopy. In general, the energy spectra of the electrons emitted 
are quite complicated because many features appear in such spectra, related 
to the different collisional processes involved before low-energy secondary- 
electron emission. As a consequence, a better understanding of the collisional 
events occuring in the surface layers before emission should allow a more 
general understanding of surface physics, including, for example, plasmon 
excitation. 

When a particle beam (with energy exceeding some threshold value) im- 
pinges on a solid target, it stimulates the emission of secondary electrons 
through collisions with target atoms. On the other hand, a fraction of the 
particles of the primary beam is also ejected from the surface because some 
particles come back and emerge from the surface after a number of elas- 
tic and/or inelastic collisions with the target atoms. If the target is not a 
thin film (i.e. if there are no transmitted particles), the remaining primary 
particles are trapped in it. The ratio between the numbers of backscattered 
and total particles (backscattered -I- trapped) is generally called, for bulk 
targets, the backscattering coefficient. If the primary particles are electrons, 
then the spectrum of the secondary electrons is clearly contaminated by the 
contribution of the backscattered primary electrons. On the other hand, as 
recently noted by Overton and Coleman [17], the problem of distinguish- 
ing between true secondary electrons and backscattered electrons is absent 
if the secondary-electron emission is stimulated by a positron beam. The au- 
thors of [17] performed an interesting experimental study of the spectra of 
fast secondary electrons uncontaminated by backscattered electrons. In fact, 
analysis was carried out for secondary electrons produced by positron beams 
impinging on a copper target at a glancing angle and at a 35° incidence 
angle, for primary energies in the range from 50 eV to 2 keV [17]. Our re- 
sults concerning the differential spectra of the energy distribution function of 
the secondary electrons emitted from a copper target above the low-energy 
cascade peak always show a clear linear trend, in agreement with both the 
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Fig. 6.9. Energy distribution of the secondary electrons emitted from Cu stimu- 
lated by a 300 eV positron beam. -I-: Monte Carlo data [18]. •: Overton and Coleman 
experimental data [17] 




E (eV) 

Fig. 6.10. Energy distribution of the secondary electrons emitted from Cu stim- 
ulated by a 1000 eV positron beam. +: Monte Carlo data [18]. •: Overton and 
Coleman experimental data [17] 
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E (eV) 

Fig. 6.11. Energy distribution of the secondary electrons emitted from Cu stim- 
ulated by a 2000 eV positron beam, o: Monte Carlo data [18]. •: Overton and 
Coleman experimental data [17] 



Sickafus law [19, 20, 21] and Overton and Coleman’s [17] experimental data 
(in the positron primary-energy range from 200 eV to 2 keV) when plotted 
on a log-log (base 10) scale [18] . The calculation of the coefficient m of 
the Sickafus law has been performed for the positron primary-energy range 
from 200 eV to 2 keV at different incidence angles without any significant 
dependence on the incidence angle being observed, in agreement with the ex- 
perimental results. In Figs. 6.9-6.11 we present our calculated energy spectra 
of the secondary electrons emitted following irradiation of copper with 300, 
1000 and 2000 eV positrons, respectively: the comparison with the experi- 
mental data of Overton and Coleman [17] demonstrates excellent agreement. 
The numerical results for m, as a function of the incident positron energy, 
together with the Overton and Coleman experimental data, are reported in 
Fig. 6.12. The incidence angle used in the calculation was 35° with respect to 
the surface of the sample. The results show excellent agreement for positron 
primary energies higher than «100 eV. However, when the positron primary 
energy is 100 eV or less, there is no agreement between theory and experi- 
ment. We attribute this discrepancy to the fact that the energy loss cannot 
be regarded as continuous in such a low-energy regime. In fact, to describe 
energy loss in such a low energy-regime, more appropriate calculations are 
necessary (quantum Monte Carlo method). 
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Eo (eV) 

Fig. 6.12. Sickafus index m for secondary-electron emission from Cu stimulated 
by positrons, as a function of the positron primary energy Eq. -f: Monte Carlo 
data [18], including numerical errors due to statistical sampling. •: Overton and 
Coleman experimental data [17] 
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A Matrices and Operators 



A.l Representation of Linear Operators 

It is always possible to represent a linear operator L as a matrix. The repre- 
sentation of the linear operator L on the basis represented by a complete set 
of eigenfunctions {rt„(r)} is given by the matrix constituted by the following 
matrix elements: 

Lnm — i'^TU LUm) — j L Um d X . (A.l) 

If L = L\ then = Lmn- In detail, 

— ^mn ' 

If {u„(r)} is an orthonormal set of eigenfunctions of the Hilbert space, then 

{Un^Um) — ^nm ■ ('^■^) 

If {rt„(r)} is a set of orthonormal eigenfunctions of the operator L, then the 
representation of L on the basis {rt„(r)} is a diagonal matrix. This can be 
written as 

LUn — ; (-^■^) 

and, as a consequence, 

d'nm — — Xn^nm • (-^■^) 



A. 2 Matrix Transformations 

Let us introduce two complete sets of eigenfunctions {M„(r)} and {u'^{r)}. 
If (fi{r) is a function in a Hilbert space, then 

‘Pir) = CnU„{r) = Y cWuir) , (A.6) 

n n 

where 

Cn — (7/7),, (A. 7 ) 
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and 

c'„ = «, V>) ■ 

We can therefore write that 

^ ^ T-fan'^m ■ 
m 

The elements of the matrix T can be calculated from 




(A.8) 

(A.9) 

(A.IO) 



The matrix T is unitary, i.e. TT^ = T^T = I where I represents the unity 
matrix: 

(TTt)„p = ^T,„„r;„= / f d\ d^r' u*^{r)up{r')6{r' - r) , (A.ll) 

n 

where we have used the closure relation 

= S{r - r') . (A.12) 

n 

Therefore, 

{TT^)rnp — j ^ ^ “ {’^rru ’^p) — ^mp ; (A. 13) 

or 

TT^ = I . (A. 14) 

In a similar way, we can show that T^T = I. 

Note that 

m 

In detail, 

c'n = «) t/’) = E V?) = E . (A.16) 

m m 

Let us indicate by A a Hermitian matrix expressed in the basis w„, and 
by A' its representation in the basis u'^^. We can write that 

A' = T^AT. (A. 17) 

In detail, 

^nm “ ^'^m) ~ ^ 'X '^jnUj , ATj-mUk) ■ (A. 18) 

jk 



y^{TjnUj,ATkmUk) = '^T*^{Uj,Auk)Tkm = '^(T'^)njAjkTkm , (A. 19) 

jk jk jk 



Now 
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and therefore 

^nm — ) nj -AjkTkm ■ (A. 20) 

jk 

Different representations of the same operators are related by unitary 
transformations. These transformations do not change the norm of vectors: 

(</>,</>) = (0,TtT<))) = (T</),r0). (A.21) 

If A is Hermitian, then A! is also Hermitian: 

{A')^ = (T^AT)^ = (T^)(tU)^ = (T^A^T») = AT = A' . (A.22) 

The Hermitian matrix A can be diagonalised by a unitary transformation, 
and the diagonal elements are matrix eigenvalues, which are real and can be 
obtained from the following equation: 

det(A - A/) = 0 . (A.23) 



A. 3 Commuting Operators 

Let us show that if two operators have a complete set of eigenfunctions in 
common, then they commute. Let Ui be a complete set of eigenfunctions that 
the operators A and B have in common. If Oj and bi are the eigenvalues, 
Aui = aiUi and Bui = biUi. As [A, B]ui = (uibi — biai)ui = 0 and 

[A,B]ip=[A,B] (A.24) 

i 

then as a consequence, 

[A,B]ip = 0 (A.25) 

for every function in the Hilbert space. This means that \A, B] = 0. 

Let us now assume that the operators A and B commute. Let be a an 
eigenvalue of A and N the dimension of the eigenspace corresponding to a, 

i.e. 

Aui = aui , (A. 26) 

for each i = 1, ... ,N. As ABui = BAui, then A{Bui) = a{Bui), so that the 
vectors Bui belong to the eigenspace corresponding to a, or 

N 

Buj = CjiUi . (A. 27) 

i=l 

When = 1, we can write 
Au\ = au\ 



(A.28) 
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and 



Bui = ciiUi , (A. 29) 

so that, as a consequence, A and B have a complete set of eigenfunctions in 
common, represented by the single eigenvector ui. 

In general, if 

Cji = {ui,Buj) (A. 30) 

and B is Hermitian, then 

c*i = (Buj,Ui) = (uj,Bui) = Cij , (A. 31) 

so that there is a unitary transformation s such that the matrix d defined 
by s^cs = d is diagonal. If dj are the diagonal elements of d, they are real 
because d^ = d. So we can write that 

N N 

^ ^ ^ j ^ (A. 32) 

i=i i=i 

and, as a consequence, 

N 

Y^s*^c,, = dkS*^. (A.33) 

i=i 

Utilising (A. 27) and (A.33), we can write 

N N N N N 

^ = XI = X X = 4 X ■ (^-34) 

j=l j=l j=l i=l i=l 

Defining 

= X , (A.35) 

we can conclude that 

Avk = avk (A. 36) 

and 



Bvk = dkVk , 

so that A and B have a basis of eigenfunctions ({'Cfc}) in common. 



(A.37) 



B The Dirac Notation 



B.l Ket and Bra Vectors 

In the Dirac notation, we use the symbol |it) to represent vectors. These are 
called “ket vectors”. Linear combinations of kets are kets. A set of kets is 
said to be linearly independent if none of the kets is a linear combination of 
the others. When a vector space has exactly N linearly independent vectors, 
it is said to have N dimensions. 

The dual space is the space of the “bra vectors” , which are indicated by the 
symbol (w|. The symbol (m| indicates the bra vector corresponding to the ket 
vector |m). The relation between a ket vector and the corresponding bra vector 
is the same as that between a wave function and its complex conjugate. The 
scalar product is represented by the bracket (|) and possesses the following 
properties: if |m) and |w) are two kets, then their scalar product is the product 
between the bra vector (u| and the ket vector jr:), and {u\v) = {v\u}*. The 
norm of |u) is real and non-negative and is given by the following: 

hf = Hu). (B.l) 

{u\u) is zero if and only if it = 0. Each linear correspondence between two 
kets I It) and |i:) defines a linear operator A such that 

li;) = A\u) . (B.2) 

Let us now consider a complete set of eigenvectors \n) such that 

\v) = '^Cn\n) = ■ (B-3) 

n n 

The coefficients c„ of the expansion are given by 

Cn = {n\v) . (B.4) 

Consequently, 

^\u) = '^\n){n\A\u) . (B.5) 

n 

The set \n) is said to be a complete set of orthonormal eigenvectors if both 
the orthonormality condition, 

(7l|in) — — ^nm ; 
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and the closure relation, 

J2\n){n\ = I, (B.7) 

n 

hold. 

Let us consider a generic vector \(p). This vector can be expanded in a 
linear combination of eigenvectors. This can be written as 

^ \n){n\ip) = I ^ |n)(n| j \ip) = I\^p) = \^p) . (B.8) 

Let us calculate the matrix elements of a linear operator: 

L = \n){n\L^\m){m\ = |n)(m|(n|L|m) . (B.9) 

n m n m 

As a consequence, we can write that 

Lnm = (n|L|m) . (B.IO) 

Let us now define the adjoint operator of L, L^, by 

{u\L\v)* = (v\L^\u) . (B.ll) 

If |w) = A\u), then (wl = {u\A^ . In detail, 

(^;| = Y,{n\A\u)*{n\ = ^(u|AV)(n| = . (B.12) 

n n 

A projector P is defined as a linear operator such that 



PP = P. (B.13) 

The projector on the direction of the normalised ket ju) is 

Pu = \u){u\. (B.14) 

If |</>) is a ket vector, then 

Pu\<i)) = |■u)(■u|(/>) = ((m|0))|'u) . (B.I5) 

Let us now show that Pu is a projector: 

PuPu = |'u)('u|m)(m| = |m)(m| = Pu ■ (B.I6) 



B.2 Continuous Spectrum 

Let us consider the eigenvalue equation that defines the spectrum of the 
position operator x, 

X(pa{x) = a(fa(x) , 



(B.17) 
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with eigenfunctions given by 

(fia(x) = S(x - a) . (B.18) 

The spectrum of the positions of a particle is characterised by a continuous 
set of eigenfunctions. They are S distributions, i.e. (intuitively) functions that 
are zero in all cases apart from when x = a, where they become oo. 

If we introduce the orthogonality relation for a continuous spectrum, the 
scalar product between two eigenfunctions belonging to two different eigen- 
values a and b is given by 



(a|6) = {(pa,‘Pb) 



dx 6{x — a)S{x — b) = S{a — h) . 



(B.19) 



The closure relation for a continuous spectrum is given by 



J da (pa{xi)(pa{x 2 ) ~ J — a)5{x2 — a) = 5(a;i — X 2 ) ■ (B.20) 



Let us consider a wave function (p{x), 

^W = /dx'i(.-xW). 

and apply the closure relation 

ip{x) = J dx' j da ^Pa{x)(fa{x')^p{x') . 

Then 

(p(x) = y Ca(fa{x) , 

where 

Ca = y dx' (pa{x')(p{x') = = {a\(f) = (p{a) . 



(B.21) 

(B.22) 

(B.23) 

(B.24) 



The coefficients of the expansion are given by the wave equation and, as a 
consequence, |(/?(a;)p is the probability density of the position operator x. Let 
us now calculate the matrix elements for the position x and the x component 
of the momentum p = {Ti / i){d / dx)\ 



Xab = (a\x\b) 



{(fiajXpb) 



J dx S{x 



a)x5{x 



b) 



aS(a — b) (B.25) 



and 



Pab 



= {a\p\b) = {(pa,PPb) = J dx S{x - a)^-^ S(x - b) 



i da 



S{a — b) . 



(B.26) 
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B.3 The Schrodinger Equation in the Dirac Notation 



The Schrodinger equation in the Dirac notation is given by the following: 

. (B.27) 

From (B.27), let us show that we can obtain the Schrodinger equation in the 
usual notation. Note that the wave function is the projection of \iIj) on (a:|: 

{x\ip,t) = . (B.28) 

From (B.27), we obtain the result that 

(B.29) 

and, taking the closure relation into account, 

= l{x\H\il),t) = J d^{x\H\^{^\'ilj,t) . (B.30) 

(ccj is independent of time, and 

{x\H\0 = {x\^ + TIC) = ~^^S{x - C) + V{x)S{x - 0 ■ (B.31) 
2m 2m ox^ 

Therefore, 









il^{x,t) 



(B.32) 



This is the Schrodinger equation in the usual notation. 



C Special Functions 



C.l Legendre Polynomials 

and Associated Legendre Functions 



The Legendre polynomial of degree I {I = 0,1; 2, 
formula 



Piiu) 



1 d‘ 






ly. 



. , oo) is defined by the 



(C.l) 



It is a polynomial with I zeros in the range (—1, +1) and with parity (— 1)^ 
The first five Legendre polynomials are 



Po = 1 , (C.2) 

Pi=u, (C.3) 

P2 = i(3«2-l), (C.4) 

Fa = i(5u^-3u) , (C.5) 

P4 = 1(35m^-30m^ + 3) . (C. 6) 

8 



Let us now introduce the associated Legendre functions: 

rjm 

pr{u) = {i-u^y^/^^^—Pi{u). (C.7) 

The Legendre polynomials are the particular associated Legendre functions 
corresponding to m = 0: 



Pi{u) = P°(w) . (C.8) 

The associated Legendre functions satisfy the orthonormality relations 



21 + 1 {I — m)l 



Ski , 



(C.9) 
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and the differential equation 



< 1 - “^ 4 -^ 4 + 



1 — 



P[^ = 0 



(C.IO) 



The following recursion relations are very useful when using the Legendre 
polynomials: 



(l-m+ l)P[!^,{u) + {l + = {21 + l)uPr{u) , (C.ll) 

(1 - u^)^Pr{u) = {l + m)Pl^,{u) - luPriu) . (C.12) 



In order to find the recursion relations corresponding to the Legendre poly- 
nomials, we can impose m = 0 in the previous equations and obtain the 
following: 

{I + l)Pi+i{u) + lPi-i{u) = {21 + l)uPi{u) , (C.13) 

(1 - u^)^Pi{u) = lPi-i{u) - luPi{u) . (C.14) 



The Legendre polynomials are eigenfunctions of the square of the orbital 
angular momentum: 



= -h 



2 



1,902 + sin2 0 902 ) 



(C.15) 



L^Pi{cos9) = Ti^l{l + l)P;(cos0) . 



(C.16) 



C.2 Bessel Functions 

Let us introduce the Bessel equation of order v. 

The solution of this equation is a linear combination of the Bessel functions 
J-i, and J+i,. 

Let us now consider the Schrodinger equation of a particle in a constant 
potential Vq, 

{V"^ + - Uo)P = 0 , (C.18) 

where V = {d/dx,d/dy,d/dz), K = p/m, = 2mE/Ti^ and Uq = 
2mVo/fi‘^. In order to proceed, let us expand the wave function in Legen- 
dre polynomials: 

E{r,cos9) = g; P; (cos 0) . 



(C.19) 
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Taking into account (C.16), the Schrodinger equation becomes the following: 
2d l{l + l) 



E 



dj.2 



■ dr 



+ - Uo 



a;E!lp,(cos6») =0 . (C.20) 



All the coefficients of the expansion must then satisfy the following differential 
equations: 



2d l{l + l) 



dr"^ r dr 
If we define 

e = K^-Uo, 

we can write 

_ /(^ + 1 ) 

dr2 J.2 



+ K^-Uo 



a,^ = 0. 






yi{r) = 0 . 



(C.21) 

(C.22) 

(C.23) 



Let us now introduce the variable x = kr, so that the last equation can be 
rewritten as 

r (P l{l + l) 



dx'^ 



+ 1 



yi{x) = 0 . 



The spherical Bessel function of order I is defined by 



ji{x) = . 



(C.24) 



(C.25) 



Since x^P is a solution of (C.24), we can conclude that the function 
krji{kr) is a solution of (C.23). 

In a similar way, once the spherical Neumann function of order I (also 
known as the irregular spherical Bessel function of order 1) has been defined 
as 



m{x) = {-iy+^^I^J_i_i/ 2 {x) , 



it is possible to show that krni{kr) is a solution of (C.23). 
The first three regular spherical Bessel functions are 



sm X 

Jo = , 

X 

sin X cos X 
Ji = — 



J2 = 



— sm X cos X , 



3 

x^ 



while the first three spherical Neumann functions are 
cosx 

no = , 



(C.26) 

(C.27) 

(C.28) 

(C.29) 

(C.30) 
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ni = -- 



cos X sin X 



1 



U2 = { T H — cos X sin X . 



It can be proved that 



3i{x) 

and that 
ni{x) 



ai— 



1 • 3 • • (21+ 1) ’ 






1C— 

The following equation holds: 

ji(0) = 5io ■ 



(C.31) 

(C.32) 

(C.33) 

(C.34) 

(C.35) 



The asymptotic behaviour of the Bessel and Neumann functions is described 
by the following equations: 



, . sin(x — Itt/2) 

JI{X) ^ 

X 

1C— ^oo 


(C.36) 


, , cos(x — ^7r/2) 

ni{x) ~ 


(C.37) 



If we indicate by fi any linear combination of the Bessel and Neumann func- 
tions (/; = aji + hni, where a and b are arbitrary coefficients), we have 



xfi-i - {21 + l)fi + xfi+i = 0 
xfi-i -{1 + l)fi - x^ = 0 . 



(C.38) 

(C.39) 



C.3 The Spherical Harmonics 

The spherical harmonics are the eigenfunctions common to the operators 
L'^ and Lz- With / = 0, 1, ... and m = —I, —l + l, ... , I — 1,1, we have 

L^Y^m ^ 

and 

LzYi"^ = hmY{^ . 



(C.41) 
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When m > 0, the spherical harmonics are given by the following: 

Yne, </>) = exp(*m^) , (C.42) 

whereas, when m < 0, the spherical harmonics can be calculated from: 

= . (C.43) 

Note that the spherical harmonics are normalised to unity on the unit 
sphere. They satisfy the orthonormality and closure relations and form a 
complete orthonormal set of square-integrable functions, as below: 

pTT p‘2tz 

/ sine d9 / d<P Y^*{9, 4>)yr{d, 4>) = SjiSkm , (C.44) 

Jo Jo 

oo I 

^ ^ Yne, <P)Yr(0', </>') = - 12') . (C.45) 

m——l 




Index 



Absorbed electrons 2, 3, 53 
Angular momentum 5, 7, 13, 100 

Backscattered electrons 1, 2, 56, 61, 
69, 87 

Backscattering coefficient 1, 53, 
55-57, 61, 65, 67, 69, 76, 84, 87 
Bessel functions 4, 27, 28, 36, 100, 101 
Bethe formula 45, 51, 54, 75 

Central potential 13, 17, 20 
Continuous spectrum 97 
Coulomb potential 20, 22 
Cross-section 3, 17, 18, 20-22, 25, 33, 
36-41, 47, 51, 54, 57, 65, 69, 70, 74, 
76, 77, 81, 82 

Density matrix 22-25 
Depth distribution 1, 3, 60, 65, 82, 86 
Dielectric function 45, 46, 76 
Dirac equation 3, 10-15, 25, 26 
Dirac notation 4, 6, 95, 98 
Dirac-Hartree-Fock-Slater atomic 
potential 4 

Eigenfunctions 91, 93, 94, 97, 100, 102 
Eigenvalues 6, 7, 12, 15, 93, 97 
Elastic mean free path 22, 74, 78 
Elastic scattering 3, 15, 17, 18, 20-22, 
25, 33, 36-41, 69, 74, 76, 77, 79, 81 
Energy and angular distribution 56, 
61, 69 

Exchange effect 37, 47 

Fermi energy 46, 48, 64, 81 
First Born approximation 18, 20, 37, 
74 

Hartree-Fock atomic potential 3 

Inelastic mean free path 4, 46, 47, 51, 
76 



Inelastic scattering 4, 79 

Legendre polynomials 4, 28, 29, 99, 
100 

Matrix transformations 91 
Monte Carlo 1-4, 65, 69, 70, 72, 73, 
76-79, 81-86, 88-90 
Multiple reflection method 61, 65-67, 
84, 85 

Operators 4-7, 9-14, 18, 19, 22, 23, 
91, 93, 95-97, 102 

Pauli matrices 8, 9, 11, 23, 32 
Phase shifts 3, 28, 35, 36 
Plasmons 2, 47, 48, 67, 69, 74, 87 
Poisson distribution 72 
Polarisation 22, 24, 33, 41 
Positrons 1, 3, 15, 47, 65, 84, 86, 87, 
89, 90 

Pseudo-random-numbers 73 

Random variables 3, 70-73 
Relativistic partial-wave expansion 
method 3, 37-41, 76 
Rutherford formula 21, 22, 37-40, 76 

Schrodinger equation 10, 98, 100, 101 
Secondary electrons 1, 2, 61-63, 69, 
80, 81, 87-89 

Simulations 1-4, 65, 72, 76, 77, 81, 
82, 86 

Solid state effect 37 
Spherical harmonics 4, 30, 102, 103 
Spin 3, 7, 8, 10-12, 15, 22-26, 31, 32 
Stopping power 4, 44-47, 50, 51, 69, 
70, 74, 76, 78, 81 
Surface films 2, 66, 67, 77, 85 

Thin films 1-4, 56-60, 65, 67, 69, 82, 
84 

Transmitted electrons 53, 69 



Maurizio Dapor: Electron-Beam Interactions with Solids, STMP 186, 105—105 (2003) 
© Springer- Verlag Berlin Heidelberg 2003 



springer Tracts in Modern Physics 



146 Low-Energy Ion Irradiation of Solid Surfaces 
By H. Gnaser 1999. 93 figs. VIII, 293 pages 

147 Dispersion, Complex Analysis and Optical Spectroscopy 

By K.-E. Peiponen, E.M. Vartiainen, and T. Asakura 1999. 46 figs. VIII, 130 pages 

148 X-Ray Scattering from Soft-Matter Thin Films 
Materials Science and Basic Research 

By M. Tolan 1999. 98 figs. IX, 197 pages 

149 High-Resolution X-Ray Scattering from Thin Films and Multilayers 
By V. Holy, U. Pietsch, and T. Baumbach 1999. 148 figs. XI, 256 pages 

150 QCDatHERA 

The Hadronic Final State in Deep Inelastic Scattering 
By M. Kuhlen 1999. 99 figs. X, 172 pages 

151 Atomic Simulation of Electrooptic and Magnetooptic Oxide Materials 
By H. Donnerberg 1999. 45 figs. VIII, 205 pages 

152 Thermocapillary Convection in Models of Crystal Growth 
By H. Kuhlmann 1999. loi figs. XVIII, 224 pages 

153 Neutral Kaons 

By R. Belusevi 1999. 67 figs. XII, 183 pages 

154 Applied RHEED 

Reflection High-Energy Electron Diffraction During Crystal Growth 
By W. Braun 1999. 150 figs. IX, 222 pages 

155 High-Temperature-Superconductor Thin Films at Microwave Frequencies 
By M. Hein 1999. 134 figs. XIV, 395 pages 

1 56 Growth Processes and Surface Phase Equilibria in Molecular Beam Epitaxy 
By N.N. Ledentsov 1999. 17 figs. VIII, 84 pages 

157 Deposition of Diamond-Like Superhard Materials 
By W. Kulisch 1999. 60 figs. X, 191 pages 

1 58 Nonlinear Optics of Random Media 

Fractal Composites and Metal-Dielectric Films 
By V.M. Shalaev 2000. 51 figs. XII, 158 pages 

1 59 Magnetic Dichroism in Core-Level Photoemission 
By K. Starke 2000. 64 figs. X, 136 pages 

160 Physics with Tau Leptons 

By A. Stahl 2000. 236 figs. VIII, 315 pages 

1 6 1 Semiclassical Theory of Mesoscopic Quantum Systems 
By K. Richter 2000. 50 figs. IX, 221 pages 

162 Electroweak Precision Tests at LEP 

By W. Hollik and G. Duckeck 2000. 60 figs. VIII, 161 pages 

163 Symmetries in Intermediate and High Energy Physics 

Ed. by A. Faessler, T.S. Kosmas, and G.K. Leontaris 2000. 96 figs. XVI, 316 pages 

164 Pattern Formation in Granular Materials 
By G.H. Ristow 2000. 83 figs. XIII, 161 pages 

165 Path Integral Quantization and Stochastic Quantization 
By M. Masujima 2000. o figs. XII, 282 pages 

166 Probing the Quantum Vacuum 

Pertubative Effective Action Approach in Quantum Electrodynamics and its Application 
By W. Dittrich and H. Gies 2000. 16 figs. XI, 241 pages 

1 67 Photoelectric Properties and Applications of Low-Mobility Semiconductors 
By R. Konenkamp 2000. 57 figs. VIII, 100 pages 




Springer Tracts in Modern Physics 



168 Deep Inelastic Positron-Proton Scattering in the High-Momentum-Transfer Regime of HERA 
By U.E Katz 2000. 96 figs, VIII, 237 pages 

169 Semiconductor Cavity Quantum Electrodynamics 

By Y. Yamamoto, T. Tassone, H. Cao 2000. 67 figs, VIII, 154 pages 

170 d-d Excitations in Transition-Metal Oxides 

A Spin-Polarized Electron Energy-Loss Spectroscopy (SPELLS) Study 
By B, Fromme 2001, 53 figs, XII, 143 pages 

1 7 1 High-Tc Superconductors for Magnet and Energy Technology 
By B, R, Lehndorff 2001, 139 figs, XII, 209 pages 

1 72 Dissipative Quantum Chaos and Decoherence 
By D, Braun 2001, 22 figs, XI, 132 pages 

173 Quantum Information 

An Introduction to Basic Theoretical Concepts and Experiments 

By G, Alber, T, Beth, M, Horodecki, P, Horodecki, R, Horodecki, M, Rotteler, H, Weinfurter, 

R, Werner, and A, Zeilinger 2001, 60 figs, XI, 216 pages 

174 Superconductor/Semiconductor Junctions 

By Thomas Schapers 2001, 91 figs, IX, 145 pages 

175 Ion-Induced Electron Emission from Crystalline Solids 
By Hiroshi Kudo 2002, 85 figs, IX, 161 pages 

176 Infrared Spectroscopy of Molecular Clusters 
An Introduction to Intermolecular Forces 

By Martina Havenith 2002, 33 figs, VllI, 120 pages 

177 Applied Asymptotic Expansions in Momenta and Masses 
By Vladimir A, Smirnov 2002, 52 figs, IX, 263 pages 

178 Capillary Surfaces 

Shape - Stability - Dynamics, in Particular Under Weightlessnes 
By Dieter Langbein 2002, 182 figs, XVIll, 364 pages 

179 Anomalous X-ray Scattering 
for Materials Characterization 
Atomic-Scale Structure Determination 

By Yoshio Waseda 2002, 132 figs, XIV, 214 pages 

1 80 Coverings of Discrete Quasiperiodic Sets 
Theory and Applications to Quasicrystals 

Edited by P, Kramer and Z, Papadopolos 2002, 128 figs,, XIV, 274 pages 

1 8 1 Emulsion Science 

Basic Principles, An Overview 

By J, Bibette, F, Leal-Calderon, V. Schmitt, and P, Poulin 2002, 50 figs,, IX, 140 pages 

1 82 Transmission Electron Microscopy of Semiconductor Nanostructures 
An Analysis of Composition and Strain State 

By A, Rosenauer 2003, 136 figs,, XII, 238 pages 

183 Transverse Patterns in Nonlinear Optical Resonators 

By K, Staliunas, V, J, Sanchez-Morcillo 2003, 132 figs,, XII, 226 pages 

184 Statistical Physics and Economics 
Concepts, Tools and Applications 

By M, Schulz 2003, 54 figs,, XII, 244 pages 

185 Electronic Defect States in Alkali Halides 
Effects of Interaction with Molecular Ions 
By V, Dierolf 2003, 80 figs,, XII, 196 pages 

186 Electron-Beam Interactions with Solids 

Application of the Monte Carlo Method to Electron Scattering Problems 
By M, Dapor 2003, 27 figs,, X, 110 pages 




Printing (Computer to Plate): Saladruck Berlin 
Binding: Sttirtz AG, Wurzburg 




